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1 Introduction 

The supersymmetry method pioneered by Efetov has proven to be a very powerful 
tool of study of the statistical properties of energy levels and eigenfunctions in disor- 
dered systems. The aim of these lectures is to present a tutorial introduction to the 
method, as well as an overview of the recent developments. One of the major achieve- 
ments of Efetov |]l| was a proof of applicability of the random matrix theory (RMT) to 
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a disordered metallic sample. More recently, the focus of the research interest shifted 
to the study of system-specific deviations from the universal (RMT) behavior. This 
will be the central topic of the present lectures. Since these deviations are determined 
by the underlying classical dynamics, this issue is closely related to the subject of 
quantum chaos. 

We begin by introducing the basic notions and ideas of the supermatrix cr-model 
method in Sec. H, considering the level statistics in the Gaussian unitary ensemble 
(GUE) as the simplest example. In Sec. |^ we turn to the problem of level correlations 
in a diffusive sample. We outline a derivation of the diffusive cr-model and discuss 
deviations of the level statistics from RMT. Section ^ is devoted to the eigenfunction 
statistics. Finally, in Sec. |^ we discuss very recent ideas of application of the method 
to chaotic ballistic systems. 

It is appropriate to give here references to a few recent review articles and books 
which are close in their topics to the present lectures. These are the Efetov's book 
the reviews by Guhr, Miiller-Groeling, and Weidenmiiller [^j and by the author [^, 
and the proceedings volume of the NATO Advanced Study Institute (in particular 
the review by Altland, Offer, and Simons P] there). Other references to the relevant 
literature will be given in appropriate places below. 

2 Introduction to the supersymmetry method and 
apphcation to RMT 

2.1 Green's function approach 

We consider the Gaussian unitary ensemble (GUE) defined Q as an ensemble of 
N X N {N oo) Hermitian matrices H = W with probability density 

7^(^) =Arexp|-yTri72| , (2.1) 

where A/" is a normalization factor. According to (|2.1| ), all the matrix elements of H 
have a Gaussian distribution with variance 

{HijH*,j,) = —5ii'5jj' ■ (2.2) 

The factor in the exponent of ( |2.1| ) allows one to keep the eigenvalues Ei of H finite 
in the limit A^ —>■ oo, with (Ef) = 1. We demonstrate below how the well-known 
results [0, for the average density of states and the two-level correlation function are 
derived within the supermatrix a-model method. 
The density of states (DOS) is defined as 

u{E) = ^TtS{E-H). (2.3) 
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The first step of tfie strategy is to express the quantity of interest in terms of the 
retarded and/or advanced Green's functions 



Gn,A = iE-H±tr])-' . (2.4) 

We have 

u{E) = -lmTT{E-H + tr])-'\^_.^o • (2.5) 

TTiV 

The second step is to write the Green's function (in a general product of 

Green's functions) as a functional integral. More precisely, in the case of RMT this 
will be an integral over a vector field with the number of components proportional 
to N; for the problem of a particle in a random potential (Sec. |^) the discrete index 
will be replaced by a continuously changing spatial coordinate, which will result in a 
functional integral. The simplest way to do it is to introduce an A^-component complex 
vector (pi, i = 1,2, . . . , N, so that 

^ "^^'^ ' /[d0*d0]exp{.E.,0:[(i5 + ^^)5..-^^*} ' 

where [d0*d0] = Hili d0*d0i. We have to put the imaginary unit i in front of the 



quadratic form in the exponent of Eq. ( p.6| ) to get a convergent integral; the convergence 
being guaranteed by ?7 > 0. 

The next step is to average over the ensemble of random matrices H. However, 
direct averaging of Eq. (|2.6|) is complicated by the fact that H enters not only the 
numerator but also the denominator. If there were no denominator, the averaging over 
H with the probability density (|2.1|) would be straightforward (a Gaussian integral 
over H). One possible way to get rid of the denominator is the replica trick first 
introduced by Edwards and Anderson in the context of the spin glass theory. The 
idea of the method is to introduce n species of the field, a = 1,2, . . . ,n. Then the 
denominator Z of Eq. ( p.6| ) is transformed to Z" and disappears in the limit n — 0. 
However, the replica trick turns out to be ill-founded for the cx-model approach 
to the problem of the level and eigenfunction statistics reviewed in these lectures. Q. 
An alternative method, which uses a combination of commuting and anticommuting 



PI variables instead of the n — ^ replica limit, was proposed by McKane and 



by Parisi and Sourlas [0]. The effective theories which are obtained in this way 



are invariant with respect to a transformation mixing commuting ("bosonic") and 
anticommuting ( "fermionic" ) degrees of freedom, which is conventionally referred to as 
supersymmetry.0 

A number of publications containing an introduction to the supersymmetry ap- 
proach are available; we mention, in addition to the Efetov's review |jl|] and his more 

^Note that if the cr-niodel is treated perturbatively (including the renormahzation group treatment 
which is a resummation of the perturbative expansion) , then the rephca trick is completely equivalent 
to the supersymmetry method. 

^Mathematicians also use the term "graded" (or, more specifically, Z2-graded) to characterize the 
arising algebraic structures. 
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recent book p[, the papers [|T4|, |T5|, |T6[. In particular, in Refs. |T^, [T^ a detailed exposi- 



tion of the calculation of the GUE level statistics with the supersymmetry is given. In 
our presentation we will concentrate on "ideological" aspects of the method and skip 
some technical details of calculations, which can be found in the literature cited above. 

2.2 Super mat hematics 

We describe now briefly the basic properties of anticommuting (Grassmannian) vari- 
ables and introduce the notions of the supermathematics which we will use. For a 
detailed exposition of the Grassmannian mathematics and of the superanalysis the 
reader is referred to the book (see also for a physicist's summary of the most 
important properties). We introduce the Grassmannian variables Xfc, ^ = 1; • • • ; 
which all anticommute to each other: 

XkXi = -XiXk, XkXi = -XiXt XkX*i = -X*iXk ■ (2-7) 
Note that Xk and xl to be considered as two independent variables. According to 



( |2.7|) , the square of a Grassmannian variable is zero. As a consequence, any function 
of Grassmanians, when expanded in a power series, may contain only the terms of the 
zeroth and the first order in each Grassmannian variable. For this reason, integration 
over Grassmannians is uniquely defined by the following rules: 

Xfcdxfc = I XfcdXfc = ; [dxk = fdxl = 0, (2.8) 



(2.9) 



^27r 

the differentials dxk, dxl anticommuting with each other and with the Grassmannian 
variables. Note that the Grassmannian integration is just a formally defined algebraic 
operation and the question "What is the domain of integration in Eqs. (|2.8| )?" does 
not make sense. 

Using the rules (|2.7|), (p.8|) , one can calculate a Gaussian integral over the Grass- 
mannians, 

f dxldxi . . . dx*NdxN exp{- xl^kiXi} = ^ ) • 

kl V ^/ 

An analogous integral over commuting variables would give det-\K/27i). This prop- 
erty of the anticommuting variables is the reason for introducing them: it will allow 
us to replace the Gaussian integral over the commuting variables in the denominator 
of Eq. ( |2.6| ) by an analogous Grassmannian integral in the numerator, thus solving the 
denominator problem! 

It is convenient to define the "complex conjugation" for the Grassmannians by the 
following rulesQ 

{Xkr = x*k; { x*kr = -Xk; ixkXiT = xlxl ■ (2.10) 

■^The notion of complex conjugation of Grassmannian variables is introduced for notational con- 
venience only; it allows one to make the treatment of Grassmannians similar to that of ordinary 
(commuting) variables and thus to introduce compact supersymmetric notations. The fact that two 
Grassmannian variables have been declared complex conjugate to each other is, however, irrelevant 
for evaluation of integrals over them, which are simply defined by the rules ( ^.8[ ) . 



5 



Furthermore, we introduce the notion of a supervector 



(2.11) 



( Si 1 

Xi 

Xn 

where Si are the commuting and Xi the anticommuting components, and n is an arbi- 
trary positive integer. A supermatrix has the structure 



a a 
I P b 



(2.12) 



where the boson-boson (bb) block a and the fermion-fermion (ff) block b are ordinary 
n X n matrices, while the boson-fermion (bf) and fermion-boson (fb) blocks a, p are 
nxn matrices with anticommuting entries. For example, for any two supervectors $, 
the tensor product $ ^E'^ is a supermatrix. The supertrace and the superdeterminant 
of the supermatrix ( 2.12| ) are defined as follows 



Str F = Tr a — Tr 6 , 

Sdet F = expStrlnF = det(a — ab^^ p)det~^b. 
Finally, hermitian conjugation of a supermatrix is defined by 



Ft 



a' p' 
-at 6t 



(2.13) 
(2.14) 



(2.15) 



where for the anticommuting (bf and fb) blocks the usual definition holds, at = (a*)^. 

It can be shown that with the above set of definitions, the usual properties of 
the vector and matrix algebra are valid for supervectors and supermatrices as well, in 
particular: 

• if F is a supermatrix and $ a supervector, then \E' = F$ is a supervector; 

• hermitian conjugation satisfies the usual properties 

F = $®^t ^ Ft = ^(g)$t, (2.16) 



($tF^)t = ^t^t$ ^ 
(Ft)t = F , etc. 

supertrace and superdeterminant of a product: 

StrFiFa = StrFgFi , 
SdetFiFs = SdetFi ■ SdetFs . 



(2.17) 
(2.18) 



(2.19) 
(2.20) 
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Gaussian integrals: 

j d<l>M<l>exp(-<l>tir$) = SdetK-^ (2.21) 

J d$M$$„$j3exp(-$tii'$) = {K~^)o,(3SdetK~^ . (2.22) 

It is assumed here that the boson-boson block /^bb of the matrix K defines a 
quadratic form with a positively defined real part, ReS^i^bbS > 0, so that the 
integral over the bosonic components S of the supervector $ converges. Integrals 



of the type (|2.22| ) with a product of a larger number of $q,'s in the preexpo- 
nential factor can be evaluated via the Wick theorem (taking into account that 
interchanging two anticommuting variables produces a minus sign). 

The notion of a supermanifold (including integration and change of coordinates 
on it) will become important for us in the course of calculation of the DOS-DOS 
correlation function (Sec. |2.4{ ). We restrict ourselves to saying that the corresponding 
definitions are natural extensions of those for ordinary analytic manifolds and refer 
the reader to Berezin's book |]T7[ (see also for a summary). A discussion of the 
structure of supermanifolds relevant to the supersymmetry treatment of the random 
matrix ensembles and a more extended list of the related mathematical literature can 



be found in 19 



Now we return to the problem of the RMT level statistics. 

2.3 Average DOS from supersymmetry 

We have, instead of Eq. ( p.6|) . 



{E - H)-,l = -I j [d$M$]5,5r exp{z ^^AES^, - H,,]^,} , (2.23) 



where 



I Xi , 



is a two-component supervector, and we have included an infinitesimally small imagi- 
nary part -|-z?7 in the definition of so that \mE > 0. The averaging over H is now 
straightforward. 



exp(z5^<|.l/f,,<|.,)\=exp|-^5^(<|.l$,)(<|.j^ . (2.24) 

\ ij / L ij ) 



The next step of our strategy is to decouple the $^ term in Eq. ( p.24| ) by introducing 
an integration over a 2 x 2 supermatrix i?, 

^^p|-^E('^l'^^)('^5'^^)} = / di?exp|-^Stri?2-z^$ti?$,| . (2.25) 
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In order that the Gaussian integral ( p. 251) over the commuting components of the 
matrix R be convergent, the latter has to be chosen in the formQ 



R 



% P 
P m ) 



(2.26) 



In other words, integration over the ff-component of the matrix R has to be performed 
along the imaginary axis. Indeed, due to the factor i in the fermion-fermion element, 
the quadratic form Stri?^ = Q'b + ?f is positively defined, ensuring the convergence. 
Substituting {^1^, in (^) , we get 



(Tt{E-H)-^) = [ [d$*d<l>] J2 SkS; exp I iE ^ <I>J<l>i -—StiR^-iJ^ [ ■ 

(2.27) 

Now we see what was the idea of the Hubbard- Stratonovich decoupling: the integral 
over the <I>-field, if taken first, is now Gaussian (and convergent due to Imii^ > 0). 
Furthermore, since the matrix R does not has any structure in the A^-dimensional 



space (where the matrices H act), the integral over each of 
produces the same factor Sdet(£' — R). As a result, we find 



1,2, 



(Tr(E - H)-^) = N j dR{E - R)^^ exp{-NS[R]} 
S[R] = -StiR^ + Strln(E - R) . 



.,N) 

(2.28) 
(2.29) 



The presence of the factor ^ 1 in the exponent of (p.28| ) allows us to apply the 
saddle-point method. The saddle point equation for the action (|2.29|) has the form 



R = {E-R)-^ 



(2.30) 



We look first for diagonal solutions of ( |2.3(J| ). According to ( |2.3(J| ), the eigenvalues r of 
the matrix R can take two values r = -^±2^/1 — E'^/A, yielding four diagonal solutions^ 



R 



E 



Si ) 



Sb, Si = ±1 



(2.31) 



Since these matrices do not belong to the original integration manifold ( p.26|) , a shift of 
the integration contours over _Rbb = % and -Rg = iq{ is necessary. However, at qy, = E 
the integrand of ( p.28| ) is singular. 



exp{-A^Strln(E - R)} = Sdet'^^E - R) oc {E - 



OO 



(2.32) 



It does not play any role whether p and p* are considered to be complex conjugate of each other 



or just two independent Grassmannian variables [see the footnote preceding Eq. (2.10)]. 

^it is implied in ( ^.3l| ) that the first term is a constant {E/2) times unit (super-)matrix. Likewise, 
we omit the unit matrix symbol in other formulas; e.g., in the l.h.s. (r.h.s.) of (2.28) E implicitly 
includes the N x N (resp. 2x2) unit matrix. 
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Recalling that ImE = +7] > 0, we conclude that integration contour over can only 
be shifted to the half-plane Imgb < 0, so that only the saddle points with Sb = +1 
are relevant. For this argumentation is not valid, since according to the definition 
of the superdeterminant the integrand has at iq^ = E a, zero [oc {E — iq^)^] rather 
than a singularity. Nevertheless, it turns out that Sf = +1 is the proper choice for the 
fermionic sector as well. The reason is, however, different: the Sf = — 1 saddle-point, 
though being legitimate, gives a contribution suppressed by as will be explained 
below. 

The leading (at ^ 1) contribution is thus given by the vicinity of the saddle 
point 

i?o = |-Vl-^V4. (2.33) 

To calculate the integral in the saddle-point approximation, we need the action at the 
saddle-point and the quadratic form around it, 

S[R,] = , (2.34) 
5''S[Ro] = C{5ql + 5q} + 2p*p) ■ C = 2-e(^- i^/l^WfA^ . (2.35) 

The Gaussian integration around Rq yields 

j dqbdqfdp*dpexp{-NC{Sql + Sqf + 2p*p)} = SdetNC = 1 . (2.36) 

[Superdeterminant of the unit matrix multiplied by a number [NC) is unity, since the 
contributions of bosons (tt/NC) and fermions {NC/tc) cancel each other.] Finally, the 
preexponential factor in (|2.28|) can be evaluated at the saddle point. 



E / E"^ 

{E - i?o)bb = (^o)bb = Y - ^1 - X • ^^-^^^ 



According to (|2.5|) , we find thus the density of states 



which is Wigner's semicircle law. 

Now we return (as was promised in the paragraph below Eq. ( |2.32| )) to another 
possible choice of sign of Sf in ( |2.31|) , i.e. Sb = — Sf = 1. The corresponding diagonal 
saddle point, 

generates in fact, by means of rotations with a Grassmannian generator, a whole man- 
ifold of saddle points. 
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The quadratic form of the action around -Ri does not contain Grassmannians, 

5^Sm = C{5ql + 5q}) , 

since the action is the same on the whole manifold. Therefore, the above compensa- 
tion of the bosonic gaussian integral (yielding a factor oc by the fermionic one 
(producing a factor oc N) does not take place, and the result is suppressed by and 
can be neglected. We will see in Sec. ^]4| that in the problem of the level correlation 
function, a manifold of saddle points emerges as well; in contrast to the present case, 
however, the balance of fermionic and bosonic degrees of freedom will be preserved, so 
that the result will be of the leading order in 



2.4 Level correlations. 

The two-level correlation function characterizing the probability density to have two 
levels separated by an energy interval u is defined as 

R2{^) ^j^^^yp • (2.39) 

As was found by Dyson and Mehta (see 0, §]), for the GUE the two- level correlation 
function has the form 

R,is) = 5{s) + l-'^^p^, (2.40) 

where s = uj/A, and A is the mean level spacing, A = 1/N {v{E)). Note that we 
assume that uj is much smaller than the energy band width; in particular, we neglect 
the change of the mean level density on the scale of uj. This is justified by the fact that, 
according to ( |2.40| ), the scale for correlations is set by the mean level spacing, A oc 
The delta function in the r.h.s of Eq. ( ^.40[ ) corresponds to a "self-correlation" of an 
energy level, the second term (unity) is a disconnected part of R2 (corresponding to 
the absence of correlations), and the last term is the non-trivial part of R2 describing 
the correlations of different levels. We show below how Eq. ( p.40| ) is obtained in the 
framework of the supersymmetric a-model method. 



We follow essentially the same strategy as was outlined in Sec. |2.3| for the aver- 
age DOS. The level corelation function is expressed in terms of the Green's functions 
as follows 

R,{uj) = ^^^^Re (TrGr-/^Tr(Gg-/^ - gI^/')) . (2.41) 

The product of two retarded Green's functions here is trivial; it can be calculated in 
the same way as the average DOS above, yielding 



-2 



(TrG^+"/'TrGg""/') ~ N-^iTiCiy ={e/2- i^l-E'^j^^ . (2.42) 



In other words, the average of the product of retarded Green's functions decouples into 
the product of the averages (the same is valid for the product of the advanced Green's 
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functions, of course). All the non-trivial information about correlations is contained 
therefore in the product of the type (GrGa), 



(2.43) 



where 



Ei^2 = E±(^^+ir]^ =E±r/2N ; Imr > . (2.44) 
To represent T2{Ei,E2) as a superintegral, we introduce supervectors of a double 



(2.45) 



size, 

r s^{{) ] 

^2(0 
I X2{i) ) 

with the first two components corresponding to the advanced and the last two to the 
retarded sector. As before, the S'-components are commuting, while the x-components 
anticommuting; the index i running from 1 to corresponds to the vector space in 
which the matrix H acts. We have then 



T2{EuE2) 



iV2 



[d$*d$] J]5*(2)^i(z)^*(j)52(j) exp{-5[$]} , (2.46) 



where the action S[^] is given by (to make notations more compact, we combine all 
S'i(z) into an A^-component vector 5*1, and the same for 5*2, Xi? X2) 



zSliE, - H)S^ - ix\{E^ - H)xi + iSl{E2 - H)S2 - ixl{E2 - H)x2 



(2.47) 



Here the matrices A and L are defined as A = diag(l, 1, —1, —1), L = diag(l, 1, —1, 1), 
with ordering of components according to ( |2.45D , i.e. (Rb, Rf, Ab, Af), where R, A 
correspond to the retarded-advanced and b, f to the boson-fermion decomposition. 
After averaging over the ensemble of matrices H, the action acquires the form 



sm = -^<^>^L [e + ^a) $ + ^ 5^(<|.jL$,)($]L<|.,) . 



(2.48) 



Decoupling of the quartic term via the Hubbard-Stratonovich transformation requires 
the introduction of a 4 x 4 supermatrix variable i?. 



T2{E,r) 



Ar2 



[d^*d^]J2smS,{t)S*2ij)S2ij) / dR 



exp |-z$tLV2 (^E-R+ ^a) _ yStri^^j 



(2.49) 



The next step of our program is to interchange the integrations / dR and J[d$*d$]. 
This leads, however, to a set of requirements of convergence of both the ^-integral 
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(when it is taken first) and the i?-integraL It turns out that in order to satisfy these 
requirements, one has to choose a non-trivial manifold of integration over the matrix 
R. Indeed, a naive attempt to generahze Eq. ( p.26| ) straightforwardly by choosing 



in boson-fermion decomposition (2.50) 



with Hermitian matrices -Rb and Rf fails, since then the ^-dependent part of the action, 
i(^'^L^/'^[E - i?)L^/^$, contains the term 

'S'i-Rb;12'S'2 + 5'2-Rb;21'S'l = 2Re 5'}i?b;12'S'2 , 

which is real and has an arbitrary sign, thus leading to a divergent integral over 5*1, 5*2. 
The natural idea to cure this problem by multiplying the components -Rb;i2 and -Rb;2i 
by i is immediately seen to fail as well, since it leads to a divergent /^-integral. The way 



out was found by Schafer and Wegner ||20[ for the case of the bosonic replica model. 



Since the problem is pertinent to the bosonic sector of the super symmetric model. 



the idea of Ref. |20| can be straightforwardly generalized to the supersymmetric case. 
For a thorough discussion of this issue, the reader is referred to the review paper by 
Verbaarschot, Weidenmiiller, and Zirnbauer 113 . 



The solution of the problem, which we are going to formulate now, is suggested by 
the invariance of the quartic term in ( [^.4(j| ) with respect to the rotations $ — >■ 
where the matrices V satisfy V'^LV = L, thus forming a pseudounitary supergroup^ 
U{1, 1|2). The reason for the above problem with convergence of the integral over the 
commuting components of $ lies in the fact that the set of hermitian matrices R^ is 
not invariant with respect to the pseudounitary group U{1, 1). The proper integration 
manifold is to be chosen as follows ^ [1^ (note that we return to the ordering 
of components according to Eq. ( p.45|) , so that the external block structure shown 
explicitly in ( ^.51[ ) corresponds to the retarded-advanced decomposition; each block 



being a 2 x 2 supermatrix) : 

i? = r[^^-^^° It-, (2.51) 

P2 + ido ) 

where 5o > is a constant (which is to be chosen below in such a way that the 
integration manifold passes through the saddle points) and Pi = P/, P2 = P^ are 
hermitian supermatrices. Further, the matrix T satisfies T'^LT = L and belongs to the 
coset spaceQ [/(l,l|2)/f/(l|l) x f/(l|l). This coset space is obtained from the group 
U{1, 1|2), if the elements of this group which commute with A (they form the subgroup 

^Let us remind that L = diag(l, 1, —1, 1); in the notation U{1, 1|2) the content of the brackets to 
the left of the vertical bar refers to the (+,— ) metric in the bosonic sector, while that to the right 
corresponds to the (+, +) metric in the fermionic sector. In other words, the group U{1, 1|2) represents 
a product [/(l,l)|bb x C^(2)|ff, "dressed" by Grassmannian generators transforming the commuting 
components into anticommuting and vice versa. 

^For a group G and its subgroup K the space G/K is formed by the set of (left) cosets gK with 
geG. 
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f/(l|l) X are identified with unity. The matrices of this coset space can be 

parametrized in the following way 



'^12 



T = 
a 

. V*2 



-irji 
ih* 



^21 



1 1/2 



t 



12 



''21 



a' 



-ri2 
ih 



\h?<l 



(2.52) 



(2.53) 



The matrix R has 16 real degrees of freedom (the same number as a 4 x 4 hermitian 
matrix); in the above parametrization 8 degrees of freedom are contained in matrices 
Pi^2 and the remaining 8 in ti2,^2i- The measure di? takes in this parametrization the 
following form (see |H] for details) 



dR = ^(Pi, P2)dPidP2d/i(T) , 



(2.54) 



where dP, are conventional matrix measures (product of differentials of all elements) 
and J-'{Pi, P2) is a function depending on eigenvalues of Pi— i^o and P2+iSo and equal to 
unity at the saddle point manifold, for which Pi 2 are proportional to the unit matrix. 
Further, d/i(T) is the invariant measure on the coset space; in the parametrization 
( P32| ), ( P35| ) of T it is given explicitly by 



d/i(T) = dti2dt2i 



(2.55) 



Having specified the manifold of matrices R over which the integration in ( p.49| ) 
goes, we can interchange the order of integrals and evaluate the integral over $ first. 



The result reads, similarly to ( p. 28 
T2{E,r)-- 
S[R] = ^StrP^ + Strln (^E-R+ 



dR E 



-1 



bb,ll 



E-R 



-1 



bb,22 



--NS[R] 



■ (2.56) 
(2.57) 



As in the case of the average DOS calculation, the large factor 3> 1 in the exponent 
allows us to use the saddle-point approximation. The saddle-point equation that is 
obtained by varying the action (p.57|) P| looks identical to Eq. (p. 301) [the difference 
being that R is now a 4 x 4 supermatrix, while it was 2 x 2 in Sec. |2.3|| . Similarly to 
Eq. (|2.31|) , diagonal solutions of this equation have the form 



R 



E 



i^l-E'^/A 



S2 



S3 



S4 



Sl, 



S4 = ±1 



(2.58) 



When fixing the signs of Sj, the same reasons as for the average DOS calculation are 
to be taken into account: 



*We neglect the term proportional to r/N when deriving the saddle-point equation, since it gives 
a correction ~ 0(1) to NS[R]. We will take this term into account below when integrating over the 
saddle-point manifold. 
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• bosonic sector: shift of the integration contour should not cross singularities. 
This requires Si = —S3 = —1; 



• fermionic sector: the leading contribution is found to be given by the saddle 
points with S2 = — S4 = —1 and S2 = — S4 = 1; contributions of the two other 
saddle-points are suppressed by 

Therefore, only 2 saddle-points out of 16 survive. In fact, it is sufficient to consider 
only one of them [say, with the signs (— , — , -|-, +)]; the second one [with (— , +, +, — )] 
will belong to the manifold obtained by rotating the first one by the matrices T. Indeed, 
if Rq = E /2 — xa/I — E'^/AA is a solution of Eq. ( p.30| ), then all matrices of the form 



R = TRqT 



E 



'I- E^/AThT-^ = 



E 



iiTuQ 



(2.59) 



are solutions [and thus saddle points of the action ( |2.57] )] as well. Here u = (y{E)) [see 
(|2.38| )] and matrices Q defined as 

Q = TAT-^ (2.60) 

(and obviously satisfying = 1) form a manifold of the supermatrix a- model. 

Comparing ( ^.59] ) with (^3T]), we see that 60 should be chosen as 60 = nu and that 
the saddle-point manifold corresponds to Pi = P2 = E/2 (and T running over the 
coset space). The Gaussian integral over Pi and P2 around this saddle-point value can 
be easily calculated, yielding unity (due to compensation of the bosonic and fermionic 
massive modes, cf. Eq. (|2.36|) ). We are thus left with an integral over the manifold 
(|2.60|) . Expanding the action ( |2.57|) up to the first order in r/A^ -C 1, we reduce it to 
the form0 

E . ^ \ f E . ^ \ ( iiTur , 



d/i(Q) 



rnvQ 



bb,ll 



bb,22 exp 



StrQA^ . (2.61) 



We have therefore reduced calculation of the level correlation function for an ensemble 
of matrices H (with ~ A^^ — > 00 degrees of freedom) to evaluation of an integral over 
a supermatrix Q parametrized by a finite number (eight) variables and belonging to 
certain non-linear space. The obtained problem is known as zero-dimensional (OD) 
cr-model. The term "zero-dimensional" distinguishes integrals of the type ( |2.61| ) over 
a single matrix Q from a field-theoretical cr-model studied in Sec. ^ with the matrix Q 
depending on spatial coordinates. 

For an explicit evaluation of integrals of the type ( p. 611) over the coset space the 
parametrization (|2.52| ), ( p.53|) is inappropriate; a much more convenient parametriza- 
tion was found by Efetov |l|]: 



Q 



Uo 



Ai 








A2 



^2 





-Ai 








-A, 



U2' 



(2.62) 



^since in Eq. ( |2.6l| ) and below we use Q as a variable on the coset space, we denote the corresponding 
invariant measure as dfi{Q); it is identical to d/i(T) introduced earlier. 



14 



Ui = exp 



' -a* 
^ a ^ 

1 < Ai < oo, -1 < A2 < 1, 



U2 = exp i 



A? 



' -(3* ] 
. P 

■1> I/^2P 



A 



(2.63) 
(2.64) 



Alternatively, instead of the variables Xi^2, /^i,2 defined in Eq. ( p.64|) , one can introduce 
the set of variables 9i, 62, 4>i, 02 via 



Ai = cosh 6^1 , fii = sinh6'ie*'^^ , A2 = cos02 , /^2 = sin 6*26*'^^ 
0<^i<oo, O<02<vr, 0< 01,2 < 27r . 



(2.65) 



Note that Ai 2 are eigenvalues of the boson-boson (and fermion-fermion) block of the 
matrix Q; we will call them simply "eigenvalues" for brevity. The integral ( p.61|) takes 
in this parametrization the form 



T2{E,r) 



E 

— - mv{\i - a*a{\i - A2)) 



X 



E 



+ z7rz/(Ai + /3*/3(Ai - A2)) 



(2.66) 



with the measure d^{Q) given by 



dAidA2 



(Ai - A2)^ 



Ld02dadaM/3d/3* 



(2.67) 



The crucial advantage of the above parametrization is that the exponent in ( p.61| ) 
acquires a very simple form (see (|2.66|) ) dependent on Ai and A2 only. This allows to 
integrate out straightforwardly the Grassmannian variables. The integration rules ( |2.8| ) 
imply that only the highest order term (~ a*a/5*/3) in the expansion of the integrand in 
a polynomial over Grassmannians gives a non-zero contribution after integration over 
da*dadp*dp. The result is easily seen to be 



2i 



dAi / dA2e^"'^'^(^i-^^) = ^sin(7rz/r)e 



f2.68) 



-1 



There exists, however, one more contribution to the integral ( |2.66| ), namely that of 
the term of the zeroth order in Grassmannians. Though naively it gives zero after 
the Grassmannian integration, the corresponding integral over Ai,2 has a singularity 
/ dAidA2/(Ai — A2)^ and thus diverges in the vicinity of Ai = A2 = 1. The reason 
for this ambiguity is in the singular character of Efetov's parametrization at Q = A, 
i.e at Ai = A2 = 1. The problem arises only for the term which does not contain 
Grassmannians; all higher order terms contain additional powers of (Ai — A2) [see 
( p.66| )], which make the integral over A's convergent. The value of the integral of the 
zeroth order term is determined by the following formula 



J d/i(g)F(Ai, A2) = F(l, 1) ^ F{Q = A) 



(2.69) 
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where F is an arbitrary function (depending on eigenvalues Ai, A2 only) which vanishes 
at Ai CX3. Equation (|2.69| ) is a particular case of a theorem (often called Parisi- 



Sourlas-Efetov-Wegner theorem in the physical literature) , which states that for a broad 
class of super symmetric models an integral of an invariant function (in our case the 
fact that F{Q) depends only on A's means that it is invariant with respect to rotations 
Q VQV^^ with V G t/(l|l) X f/(l|l)) is given by its value at the origin (in our 



case Q = K). For a discussion and a proof, see Refs. pi] , P2l . A general mathematical 
treatment of such anomalous contributions to integrals over supermanifolds can be 
found in |23|, HI. 



Applying ( |2.69| ) to the zeroth order (in Grassmannians) term of ( |2.66| ) and adding 



the result to the contribution (|2.68 ) of the highest order term, we finally get 



T2{E,r) = { ^-inu] ( ^ + inu] + sm{7rur)e''"''' . (2.70) 



2 / V 2 r 



Substituting this in ( p.41| ), including the (GrGr) contribution ( p.42|) , and taking into 



account that ur = vujN = a;/A = s, we find the two-level correlation function 

2" 



T2{E,r) - -inu 



^, , sin^fvrs) 



in agreement with ( p.40| ). 

2.5 Comments and generalizations 
2.5.1 Structure of the saddle-point manifold. 

We first note that the topology of the cr-model manifold was crucially important for 
the above calculation. Indeed, if we keep only the ordinary part (which does not con- 
tain Grassmannians) of each element Q of this supermanifold|3, we get a conventional 
(not super-) manifolcQ which is a product of the hyperboloid ?7(1, l)/f/(l) x U{1) 
(parametrized by Ai = cosh 6*1 and 0i) and the sphere U{2)/U{1) x U{1) (parametrized 
by A2 = cos 6*2 and ^2). This combination of the compact (A2) and non-compact (Ai) 
degrees of freedom is clearly reflected in the result, see Eq. ( p.68[ ). 



It is appropriate here to add a comment concerning a somewhat subtle point of 
the calculation. The non-compact (hyperbolic) symmetry of the bb-sector originates 
from the opposite signs of the first and the third terms in Eq. ( p. 47] ). This sign choice 
was unambiguously dictated by the requirement of convergence of the integral over the 
bosonic components S'1^2 of the supervector $. On the other hand, the situation seems 
to be different for the two terms of (|2.47 ) containing the Grassmannians Xi.25 since 



the integral over the anticommuting variables is always defined and no condition of 



^"^the discarded part containing the Grassmannians is called nilpotent 
called the base of the supermanifold 
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convergence arises. Therefore, the choice of signs of the second and the fourth terms in 
( |2.47|) , which has eventually led us to the compact symmetry of the ff-sector, seems to 
be arbitrary. A thorough analysis shows, however, that this apparent freedom is 
spurious, and the symmetry of the ff-block is uniquely fixed to be compact on a later 
stage of the calculation by the requirement of convergence of the integral over the coset 
space (i.e. over the matrix T). 



2.5.2 Gaussian ensembles of different symmetry. 

Let us recall that for the sake of technical simplicity we considered the case of GUE 
(with the matrix H being hermitian without any further restrictions) while calculating 
the level statistics. An analogous (though technically more involved) calculation can 
be performed for the two other Gaussian ensembles of Wigner and Dyson, the Gaussian 
Orthogonal Ensemble (GOE) of real symmetric matrices and the Gaussian Symplectic 
Ensemble (GSE) of real quaternionic matrices [0, If the matrix H is considered 
as a Hamiltonian, the GUE describes systems without the (antiunitary) time reversal 
symmetry, while GOE and GSE correspond to systems with the time reversal invariance 
realized by an antiunitary operator T such that = 1 (GOE) or = — 1 (GSE). 

Calculation of the two-level correlation function for GOE within the supersymmetry 
approach is performed essentially in the same way as for GUE, see [|T], |^ . However, due 
to the presence of an additional symmetry in the problem, one has to double the size 
of the supervectors $4, combining $j with its time reversal $*. Correspondingly, one 
ends up with a OD cr-model of 8 x 8 supermatrices Q parametrized by 16 independent 
variables. Now, instead of two for GUE, there are three eigenvalues, 1 < Ai,A2 < C)0 
and — 1 < A < 1. A similar structure is obtained for GSE with 8x8 Q-matrix as well, 
but with two "compact" and one "non-compact" eigenvalues, — l<Ai<l,0<A2<l, 
1 < A < cxD. The evaluation of the corresponding integrals |]1], reproduces again the 
result for the two-level correlation function obtained by Mehta and Dyson within the 
orthogonal polynomial methodQ 



i?2(s) = 1 + S{S) 



Sm TTS 



[ITS 



-sign(s) 



SifTTsl 



COS vrs 



sm TTS 

"MP 



R2{S) 

Si(x) -- 



l + 5{s) - 
siny 



sin^(27rs) 
(27rs)2 



+ Si(27rs) 



cos27rs sin27rs 



2713 



(27rs)2 



(GOE), 
(2.72) 

(GSE), (2.73) 



y 



dy 



In the remaining part of the lecture course, we will continue considering explicitly the 
unitary class only; corresponding results for the other two classes will be sometimes 
quoted in the end of the calculation. 



^^Note that in the symplectic case all the levels are double degenerate (Kramers degeneracy). This 
degeneracy is disregarded in ( ^.73 ), which thus represents the correlation function of distinct levels 



only, normalized to the corresponding level spacing. 
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On top of the three "standard" Wigner-Dyson ensembles, seven other ensembles 
have been introduced more recently. They become relevant, in various physical situ- 
ations, near a special point of the spectrum where the symmetry of the problem gets 
enlarged. Specifically, this happens (i) at the center of the band for a particle moving 
on a bipartite (AB) lattice with only off-diagonal (A— »-B or B— >^A) hopping allowed 
(ii) in the chiral random matrix ensembles which model the massless Dirac 



operator in the quantum chromodynamics, near the zero point of the spectrum; (iii) 
near the Fermi energy in models of a mesoscopic metallic grain in contact with a su- 
perconductor ||2^. We will not consider such ensembles in this lecture course. An 



interested reader is referred to Ref. |jT9[ for a review of related algebraic structures. 

A natural question that may be asked at this point is why is there a need in 
the supersymmetry formalism if the results discussed above have been obtained much 
earlier by the methods of "classical" RMT. In fact. Sec. will give an answer to 
this question, since the supersymmetry approach will be used there to study statistical 
properties for a much broader class of problems (disordered and chaotic systems), 
essentially different in their formulation from the RMT ensembles. However, already 
here we want to give two examples of the random matrix ensembles which, while being 
direct generalizations of the Gaussian ensembles, are not accessible within the classical 
RMT methods. 

2.5.3 Ensembles with non-Gaussian distributions of matrix elements. 

Let us consider an ensemble of large (A^ —>■ oo) N x N hermitian matrices with all 
matrix elements being independent and equally distributed, but (in contrast to GUE) 
with some non-Gaussian distribution function f{z,z*) = /(|zp). In other words, the 
overall probability density for the matrix H is assumed to be 

nH) = ll f(\H^^\") ■ (2-74) 

For the Gaussian distribution, /(|-2^|) = exp(— A^|2;p), this is reduced^^ to the Gaussian 
ensemble ( p.l| ). However, for any other distribution function / the probability density 
V{H) cannot be written in the form expTrF(ff) any more. In other words, the en- 
semble ( p^.74D is no more invariant with respect to the unitary rotations, H UHU~^ 
and the probability distributions of eigenvalues and eigenvectors of H do not decouple 
any more. This precludes the application of the orthogonal polynomial method, which 
uses an explicit form of the distribution of the eigenvalues. 



On the other hand, the supersymmetry method can still be successfully used [ p7[] . 
After averaging over H one gets, instead of the quartic term in the action ( |2.48|) , 



ij I i<j 



2\ 



^•^The distribution of the diagonal elements Ha is not important in the leading order in 
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= n {l - i^lL^J){^]L^^){\z'\)f + \{^lL^,n^]mn\z'\)f + ...[, (2.75) 

where (•••)/ denotes averaging over /(|zp). As in the case of GUE, we will assume the 
normalization which keeps the band of eigenvalues of H finite in the limit — oo. 
This implies ~ i.e. {\z'^\)f ~ If / is a smooth, non-singular function, 

then (l-z^l)/ ~ etc. For this reason, one can drop all terms beyond the quartic 

one in the second line of ( p.75| ). Indeed, the integral over $ (taken after the Hubbard- 
Stratonovich decoupling) will be dominated by $ ~ A^*^, see Eq. ( p.49| ). Therefore, the 



second line of ( |2.75| ) is in fact an expansion in 1/A^. Dropping all higher order terms 



and reexponentiating the quartic term (which is ~ 1/A^), "we get 

exp{-J2^lH,,L^,}\ ^ ]JexpL^{\z'\)fi^lL^,)i^]m)\ , (2.76) 

\ ij I i,j 

which is the same result as for the Gaussian distribution f{z, z*) = exp{ — 

Therefore, the spectral statistics for the ensemble (|2.74| ) has the same Wigner-Dyson 

form as for the GUE. 

Sparse random matrix ensemble. A notable exception from this statement is 
formed by an ensemble of sparse random matrices ^7\, |2^ for which the distribution 
function /(I^P) is singular at z = 0. More specifically. 



fi\zn = [l-j^)6{\zn + j^h{\z\'), (2.77) 

where p > is a constant of order unity (i.e. independent of A^) and h is a smooth 
(in particular, non-singular at 2; = 0) distribution function with ~ 1. The 

distribution ( |2.77|) implies that almost all matrix elements of H are zero; only p (in 
average) out of A^ elements in each row are non-zero (and distributed according to 
/i(|zp)). Because of the singular character of /(|-2p) truncation of the series in ( p.75| ) 
is no more legitimate. Indeed, we have now ~ 1/^; (k'^l)/ ~ 1/^; ^^c., so that 

all the terms have to be taken into account. To decouple all of them via a kind of the 
Hubbard-Stratonovich transformation, one has to introduce, in addition to the usual 



matrix Qa/s, also higher order tensors Q^j^^, Qal-ySfiu^ ■ ■ ■■ turns out ^ that 



they all can be combined in a function y($, $"''), so that the Hubbard-Stratonovich 
transformation acquires a functional form. The analysis show that the model exhibits 
an Anderson localization transition if one of the parameters {e.g. p) is changed. At 
p < Pc all eigenvectors are localized and the level statistics is Poisonnian (uncorrelated 
levels) in the limit N 00; at p > pc the states are delocalized and the Wigner-Dyson 
statistics applies]^; finally p = Pc is the critical point of the Anderson transition. We 
mention also that this model has effectively infinite-dimensional character and is closely 
related to the problem of Anderson localization on the tree-like Bethe lattice, see . 

^''To be precise, one should speak about the states belonging to the infinite cluster only, see [p7| [2q]. 
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2.5.4 Random banded matrices. 



The random banded matrix (RBM) ensemble is defined in the following way (see e.g. 
PU| , PT[1). As in the Gaussian ensemble, all the matrix elements are supposed to be 
independent and have a Gaussian distribution. The difference is that the variance of 
this distribution is now not the same for all the elements but rather depends on the 
distance from the main diagonal, 

{\H,,\') = F{\z - j\) . (2.78) 

The function F{r) is supposed to be roughly constant for r < b and negligibly small 
for r ^ b. Therefore, all the elements of the matrix which are essentially non-zero are 
located within a band of width ~ b around the main diagonal (hence the term "banded 
matrices"). Again, the RBM ensemble is not rotationally- invariant and thus cannot be 
treated by the classical RMT methods. In contrast, the supersymmetry approach can 
be applied to this ensemble In particular, one finds that if b/N is kept constant 



while the limit — > cxd is considered, the level statistics acquire the universal Wigner- 
Dyson form. More generally, this ensemble describes a quasi-one-dimensional system, 
with statistical properties depending on a value of the ratio X = N/b"^. If X <^ 1 (as 
in the above limiting procedure), the level (and eigenfunction) statistics are essentially 
of the GUE form; in the opposite limit X ^ 1 the system is in the strong localization 
regime. This will be discussed in more detail in Sec. |], |. 

2.5.5 Parametric level statistics. 

We have just demonstrated that the supersymmetry approach is extremely useful when 
the matrix ensemble different from the Gaussian ones are considered. However, also for 
the Gaussian ensembles the supersymmetry method has produced a bulk of new results 
for quantities more complicated than the conventional two-level correlation function 
( |2.39| ). As an important example, we mention here the parametric level statistics which 
is defined in the following way. Let us perturb the set of matrices Hq forming the GUE 
by adding some given traceless matrix V multiplied by a parameter a, 

H{a) = Ho + aV . (2.79) 

Consider now the following correlation function 

R2{uj, 6a) = (Tr 6{E - H{a)) Tr 6{E + uj - H{a + 5a))) / (yf . (2.80) 

For (5q; = this is just the usual two-level correlation function ( p.39| ). The parametric 
DOS-DOS correlation function was calculated via the supersymmetry method by Si- 
mons and Altshuler [33|. As for the simple level correlation function ( p.3£ ), the result 



has a universal form if the perturbation parameter 5a is properly rescaled, 

ms(Ai-A2)-yx2(A?-A^) , (2.81) 



R2{'^,5a) = 1 + -Re j dAi j dA2exp 



where x = 5a{{daEi)'^y^'^ / A and, as before, s = cu/A. The corresponding formulas for 
the other two ensembles can be found in pSj ^ |]. 
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3 Level statistics in a disordered sample: Diffusive 
cr-model. 



3.1 Derivation of the diffusive cr-model. 



Having introduced the main ideas and ingredients of the supersymmetry formahsm 
for the Gaussian Ensemble, we are prepared to proceed with consideration of a much 
richer problem of a particle moving in a random potential. The Hamiltonian has now 
the form|^ ^ 

H=^ + U{r), p = -zV (3.1) 

where U{r) is a disorder potential. For simplicity we choose it to be of the white-noise 
type (this is not essential for the derivation of the cr-model) 

(t/(r)f/(r')) = -^5(r-r'), (3.2) 

where r is the mean free time. We will assume that the time reversal invariance is 
broken (say, by weak homogeneous or random magnetic field) and restrict our consid- 
eration to the (technically simpler) case of the unitary symmetry.p^ The DOS and the 
Green's functions are given by a direct generalization of Eqs. (|2.3| ) - (|2.5|) , 



^^^^ = 2^/ [G'f(r,r)-Gi(r,r)] 



(3.3) 



Gg,A(ri,r2) = {v,\iE-H±tr^)-'\r2) , ^ +0 . (3.4) 

In the two-level correlation function (|2.39|) the (GrGr) and (GaCa) terms are again 
trivial (decouple in the product of averages), and all the non-trivial information is 
contained in the (GrGa) terms. Introducing the supervector field 

' Si{r) ' 
Xi(r) 

. X2(r) ^ 

we write, in analogy with (|2.46|) , (|2.47| ), the corresponding product of the Green's 
functions functional integral. 



(3.5) 



Gf+"/'(ri,ri)G 



A (r2,r2j 



D$D$t5i(ri)5*(ri)52(r2)52*(r2) 



(3.6) 



15 We set ;i= 1. 

i^We will use a somewhat sloppy (but convenient) terminology and refer to systems having the 
symmetry of GOE, GUE, and GSE as to systems of orthogonal, unitary, and symplectic symmetry, 
respectively. For a disordered electronic problem, the orthogonal symmetry corresponds to purely 
potential scattering in the absence of magnetic fields (both time reversal and spin-rotation invariance 
preserved) , the symplectic symmetry to the spin-orbit scattering (spin-rotation invariance broken, but 
time reversal preserved), while the unitary symmetry class is realized in the presence of homogeneous 
or random magnetic field breaking the time reversal invariance [^5[ 0, ^ . 
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Averaging over the disorder with the correlator (p.2|) produces a quartic term, 
exp{i j d'^r<l>t(r)Lf/(r)$(r)}^ =exp|-^-^ j d'^r [$"^(r)L$(r)]^ 



(3.7) 



the Hubbard- Stratonovich decouphng of which requires the introduction of an r-dependent 
supermatrix field -R(r), 



exp 



DR exp{-z J d'^r$t(r)Li/2^(r)Li/2^(r) 
- TTZ/r / d'^rStri?2(r)} . (3.8) 



For any given r the matrix R{r) has the structure specified by Eqs. ( ^.51| )-( [2.53| ). 
Substituting ( p.8|) , ( p.7|) in (|3.6| ) and performing the ^-integral, we obtain, similarly to 
(p.57|) , an integral over the i?-field with the action 



UJ 



S[R] = HUT / d'^r Stri?" + Str In E + -A - — - 
The corresponding saddle-point equation reads (for 0) 



(3.9) 



R{v) 



2'KVT 



9 



E-^-R 
2m 



-1 



(3.10) 



We look first for a diagonal, r-independent solution of ( |3.10| ), R = diag(gi, g2, 5'3, 94) 
which has, in the weak-disorder regime Et ^ 1, the form 



with 



^3 



(GR(r,r)) 



^Re(G.(r,r))±^, 



d'^p 

(27r)° 



(E - p'^/2m + i/2r)-\ 



(3.11) 



(3.12) 



The first term in (|3.11|) gives a non-interesting constant real contribution to R, which 
can be absorbed^] into the energy E. 

The choice of signs in the second term of (|3.11|) is dictated by the same considera- 
tions as in the GUE case (see the text below Eq. (|2.58|) ), leaving us with 



Ro = -^Re(GR(r, r)) - ^A = a - ^A. 
27rz/r 2r 2r 



(3.13) 



The manifold of saddle-points is generated from Rq by rotations with matrices T defined 
in (Pi, (Pi, 



R = a TAT-^ = a 

2r 



2r 



Q 



(3.14) 



^^As it stands, Eq. (3.12) has an ultraviolet divergence at p ^ cx) in d > 2. This is related to the 
white-noise (zero correlation length) character of the disorder potential. Assuming a finite correlation 
length of the disorder would introduce an UV-cutoff and make the integral (B.12) finite. In a realistic 
model of a disordered metal the correlation length is set by the screening length, which is of order of 
the Fermi wave length. Introducing a cutoff ^ pp = {2mE)^/'^ in ( 3.12 ), we find that the first term 
in ( 3.12 ) is much smaller than E in the considered weak-disorder regime. 
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with Q from Eq. ( p.60|) . The set of matrices ( p.l4|) constitutes a manifold of degenerate 
(at cij — > 0) saddle points of the action (|3l9|). 

We allow now the matrix T (and consequently Q) to vary slowly in space, 



R{v) 



a 



g(r) = T{T)AT-\r 



(3.15) 



While (|3.15| ) with an r-dependent Q is not an exact saddle-point any more, the fields of 
this form constitute the low-lying excitations for the functional integral j DR . . . e~^^^^ 
with the action ( p.9| ) . Performing the gradient expansion of the second term in Eq. (|3.9| ) 
for the configurations defined by Eq. ( p.l5|) (and also expanding in u up to the linear 
term), we find the action for these low- lying modes 



S[Q] 



T 



J d'^r Str[-D(Vg)^ - 2tujAQ] 



(3.16) 



where D = Vpr/d is the diffusion constant (here vp = (2i?/m)^/^ is the particle velocity 
and d the spatial dimensionality). The two-level correlation function is thus reduced 
to a functional integral over the Q-matrix field slowly varying on the coset space. 



Dg(r) 



d'^r StiQAk 



-S[Q] 



(3.17) 



where k = diag(l, —1, 1, —1), i.e k is equal to 1 (—1) in the boson-boson (resp. fermion- 
fermion) block. 

The field theory characterized by the action (|3.16| ) is called a (d-dimensional) non- 
linear (T-model. It characterizes the low-frequency long-wavelength physics of the origi- 
nal electronic problem. The fact that such a matrix a-model in the replica limit n —>■ 
is the effective field theory for the problem of a particle in a random potential was first 
realized by Wegner ; a derivation (in the replica formulation) was given by Schafer 
and Wegner |^ and by Efetov, Larkin, and Khmelnitskii The supersymmetric 
version of the model was presented by Efetov [|l|. To elucidate the physical content 
of the cr-model ( [j.l6| ), it is instructive to draw a parallel with a model of classical 
magnetic moments with a ferromagnetic interaction in an external magnetic field. 



i/[S] = -5^Jrr'S(r)S(r')-$^BS(r) 



(3.18) 



where S(r) are n-component vectors with S^(r) = 1. At low temperatures, only the 
low-energy sector is important, which corresponds to a slowly varying vector S(r). The 
Hamiltonian ( |3.18| ) is then reduced to a vector a-model. 



H[S] 



d'^r 



L2 



(VS(r))2-BS(r) 



(3.19) 



where k is the spin stiffness. Let us now demonstrate the analogies between ( ^.16| ) 
and ( ^.191 ). The ra-component unit vector S sweeps the {n — l)-dimensional sphere 
5""^^ isomorphic to the coset space 0{n)/0{n — 1). If the external magnetic field B 
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is directed, say, along ei = (1, 0, . . . , 0), then ei plays in (|3.19|) the same role as A in 
(|3.16|) . Indeed, 0(n — 1) is the subgroup of 0(n) which does not rotate ei; for this 
reason the space of matrices from 0{n) which generate different vectors when acting 
on ei is the coset space 0{n)/0{n — 1). Fully analogously, x is the 

subgroup of those matrices from f/(l, 1|2) which do not rotate (i.e commute with) A, 
so that the manifold of matrices Q is generated by rotations T belonging to the coset 
space f/(l, l|2)/f/(l|l) x [/(1|1). Furthermore, while the first term of (|3.19|) is invariant 



with respect to a global rotation of S(r), the second term breaks this 0{n)/0{n — 1) 
symmetry. In the ferromagnetic phase the symmetry is broken spontaneously at -B ^ 0; 
the corresponding Goldstone modes are the spin waves. To write their Hamiltonian 
explicitly, one should represent H in terms of independent degrees of freedom, i.e. in 
terms of the transverse (with respect to ei) part (the longitudinal component being 
^i = (l-Si)V2), 

H[S±] = ^ / h[VSx(r)]2 + 5Si(r) + 0(Si(r))] , (3.20) 

with the last term describing interaction of the spin waves. Equation ( p.20|) implies 
that the correlation function of the transverse magnetization has the Goldstone-type 
low-momentum behavior, 



/ d'^re-'^'-(Sx(0)Sx(r)) oc . (3.21) 



Comparing ( p.l9|) with (|3.16| ), we see that the frequency u plays, in the case of a dis- 



ordered electronic system, the same role of a symmetry breaking field as the magnetic 
field B for a ferro magnet. The soft modes for ( p.l6[ ) are the diffusion modes. To write 



their action explicitly, one should choose some parametrization of Q in terms of inde- 
pendent degrees of freedom. Two parametrizations, most often used for perturbative 
calculations, are 

Q = (l-H^/2)A(l-H^/2)-^ = A(1 + W^/2) (1-1^/2)-^ 



and 



A(l + W^ + — + — + — + ...) (3.22) 



Q = A(^Vr+ym^j = Ah-Fiy + ^ 



(3.23) 



with W having the off-diagonal block structure in the RA space. 
In either case the action ( |3.16| ) takes the form 



S^.|^^/d^rSt..|D(VVVT-..«- + 0(«-)]. (3,25) 
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with the last term on the r.h.s. describing the interaction of the diffusion modes. The 
analog of ( p.21| ) is the diffusion propagator (see below), 

/ d<^re-''^'-(StrA;W^i2(0)A;iy2i(r)) — ^ • (3.26) 

After this digression we return to the analysis of the two-level correlation function 



3.2 Reduction to the OD cr-model: Universal limit. 



According to (|3.25|) , the "energies" (eigenvalues of the quadratic form of the action) of 



the diffusion modes are given by 

Dq^ — iuj = — ioj , (3.27) 

with being the allowed values of wave vectors. One can show that in an isolated 
sample the boundary condition for the Q-field isQ 

V n^; I boundary 

, (3.28) 

where Vn is the normal derivative. If we consider a rectangular sample of a size 
Li X . . . X Ld, the wave vectors q^ are quantized according to 

qM = ^(^,---,^) ' ^. = 0,1,2,.... (3.29) 

For periodic boundary conditions, which are often used by theoreticians, we have in- 
stead 

qM = 27r(^^,...,^^ ' ^^ = 0,±1,±2,.... (3.30) 

In either case, one can distinguish the zero mode with qo = (and thus cq = 0) 
from all other modes with > D^n/Ly, where L = max{Li, . . . The energy 

Ec = D / is called the Thouless energy; it is of the order of the inverse time of diffusion 
through the sample. Equation ( |3.27 ) suggests that for uj <^ Ec one can neglect all non- 



zero modes when calculating the level correlation function. This approximation used 
by Efetov is known as the zero-mode approximation. As a result, the functional 
integral ( |3.17| ) reduces to an integral over a single supermatrix Q, acquiring the form 
of the correlation function of the OD cr-model, Eq. ( p.61| ). Therefore, in the zero- 
mode approximation the level statistics is described by the OD a-model and thus has 
(according to Sec. U) the RMT form (|]7l|). 

For systems with preserved time-reversal invariance the same consideration leads 
to the (T-models of the orthogonal or symplectic symmetry (depending on the 



18n 



To understand the physical meaning of ( 3.28D , one should remember that Q was introduced as 



a field conjugate to the product ^^^^ i.e. it is a kind of a density field. Equation ( 3.28 ) is thus the 
condition of the absence of current through the sample boundary, which is precisely what should be 
imposed on a boundary with vacuum or with an insulator. 
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presence of spin-dependent scattering). In the zero-mode approximation they reduce 
to corresponding OD cr-models yielding the level statistics ( p.72| ), (|2.73|) of Gaussian 
ensembles of the corresponding symmetries. 

In order for the zero-mode approximation to make sense, the Thouless energy 
should be much larger than the mean level spacing A. It is easy to see that the 
ratio g = 27rEc/A is the dimensionless (measured in units of e'^/h) conductance of the 
sample (for a current flowing in the direction of L). Therefore, the above condition 
is equivalent to (7 ^ 1; we will term the corresponding situation "metallic regime" 
(or "weak localization regime"). The opposite case, g ^ 1, corresponds to the strong 
localization of a particle and is realized in 2D or 3D if the disorder is not weak[^ 
Et ~ 1, or in sufficiently long systems of quasi-lD geometry (see below). The notion of 
the dimensionless conductance plays a central role in the scaling theory of localization, 
p9| ] for review. 



see 



3.3 Deviations from universality 
3.3.1 Perturbation theory. 

Deviations from the universal (RMT) behavior are due to the non-zero diffusion modes. 
The first calculation of the level statistics in a diffusive grain taking into account 



the non-zero modes was performed by Altshuler and Shklovskii They used the 

perturbation theory, which amounts to keeping only the terms of the leading order in 
W in the expansion of the preexponential factor and of the action ( |3.25| ) in Eq. ( |3.17|) . 
This yields (for the unitary symmetry) 

_ _ 2 

^2,As(t^)-l - jT^Re[BW{r) f d'^rStiWk 



X exp I ^ j d'^r StTW[DV^ + iu]W 
1 

= — ReV- TTT. (3.31) 

27r2 ^ (e^ - iuy ^ ' 

The same result is obtained for the other symmetry classes, with an additional factor 
2/(3, where (3 is the commonly used parameter equal to 1, 2, and 4 for the orthogonal, 
unitary, and symplectic symmetry, respectively. 

We analyze first the perturbative expression ( p.31| ) in the limit uj ^ Ec (or, equiv- 



alently, s = uiA <^ g), when the zero- mode approximation applies. Keeping only the 
term with cq = in (|3.31| ) yields 

i?2,As(^) - 1 = -77^ , s<^g. (3.32) 



To compare this perturbative formula with the exact result ( 2.71 ), we note that the 



connected part of the RMT correlator ( p.71| ) can be decomposed into the smooth and 

Strictly speaking, in 2D even for a weak disorder all states are localized, but the corresponding 
localization length is exponentially large. 
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the oscillatory contributions, 

sin^(7rs) 1 cos(27rs) 

We see that the Altshuler-Shklovskii calculation reproduces the smooth term, but fails 
to give the oscillatory contribution. This is because their method is perturbative in 
1/s, while the oscillations oc cos(27rs) are of non-perturbative (in 1/s) character. Fur- 
thermore, the Altshuler-Shklovskii result gives no information about the actual small-s 
behavior of R2{uj). 

We consider now the opposite regime, u ^ Ec {s ^ g) , in which the summation in 
(|3.31|) can be replaced by integration, yielding^ 

R2,As{oo) - 1 ^ cag-'^/h''/^-^ , s > , (3.34) 

with a numerical coefficient q. Therefore, in the domain s ^ g the level statistics 
differs completely from its universal (RMT) form. The slower (compared to RMT) 
decay of DOS correlations at s ^ g corresponds to the diffusive motion of the particle 
at times t -C t^) ~ E^^- In contrast, at long times, t 3> tn, the trajectory covers 
ergodically the whole sample volume, and the correlations acquire their universal form. 

It should be mentioned here that Altshuler and Shklovskii used in fact not the 
(T-model formalism but the impurity diagram technique, with the ladder diagrams 
representing the diffusion modes. This diffuson-cooperon diagrammatics (completely 
equivalent to the perturbative expansion of the a-model) is by now a standard tool, 
which has allowed to discover and to study a number of remarkable effects in mesoscopic 
physics, in particular the weak localization, the Altshuler-Aronov effect of interplay of 
the Coulomb interaction and disorder, and the universal conductance fluctuations. 
Several excellent reviews of these subjects are available, see Refs. [^, ^ ^ 



However, as the above discussion shows, the perturbative methods are not sufficient 
for studying the statistics of energy levels (and also of eigenfunctions, see Sec. |), and 
the non-perturbative cx-model approach has to be used. 

3.3.2 Deviations from universality at to <^ Ec. 

We return now to the region s -C (? and consider deviations from the universal (RMT) 
behavior. The method that allows us to calculate such deviations from the universality 



was developed in |44] and can be outlined as follows. We first decompose Q (in a proper 
way taking into account the non-linear constraint = 1) into the constant part Qo 
and the contribution Q of higher modes with non-zero momenta. Then we integrate 
out all non-zero modes. This can be done perturbatively provided the dimensionless 
conductance g ^ 1. As a result, we get an integral over the matrix Qq only, which has 
to be calculated non-perturbatively. 



■"We assume the spatial dimensionality d < A (otherwise, the sum in ( ^.31 ) has an ultraviolet 



divergence). We also note that in 2D the Altshuler-Shklovskii calculation gives the coefficient C2 = 0; 
a more careful consideration taking into account higher order terms of the perturbative expansion 
yields C2 ^ for the unitary and C2 ~ I/5 for the orthogonal and symplectic symmetry. 
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We begin by presenting the correlator R2{uj) in the form 
R.iuj) = (^1^ / BQeM-S[Q]}\u=o , 

S[Q] = J Str(VQ)^ + s J StrAg + uj StrQAk , (3.35) 

where s = Tis/2iV ^ u = 7Tu/2iV. Now we decompose Q in the following way: 

g(r) = ToQiv)T,-' , (3.36) 

where Tq is a spatially uniform matrix from the coset space and Q describes all modes 
with non-zero momenta. When uj -C E^, the matrix Q fluctuates only weakly near the 
origin A of the coset space. In the leading order, Q = A, thus reducing (|3.35|) to the OD 
cT-model, which leads to the Wigner-Dyson result (|2.71|) . To find the corrections, we 
should expand the matrix Q around the origin A. We use the parametrization ( p.22| ). 



Q = A(l + 1^/2)(1 - W/2)-^ = Ail + W + — + — + ...] . 



(3.37) 



To keep the number of degrees of freedom unchanged, we should exclude the zero mode 
from Q, which is achieved by the constraint J d'^rW = 0. Substituting the expansion 
(IOTI) into Eq. ( ^351) , we get 

S = So[Qo] + Sw[W] + Si[Qo,W] ; 

50 = j Str [SQoA + uQ^Ak] , 

S^ = ^j Str(Viy)^ , 

51 = ^J Str[sf/oAiy2 ^ uUokAW^] + 0{W^) , (3.38) 

where Qq = TqATq^, Uq = Tq^ATq, U^k = T^^AkT^. Let us define S'efjfQo] as a result 
of elimination of the fast (non-zero) modes: 

where (. . .)w denote the integration over W with the Gaussian weight exp{— S'h^[H^]}, 
and is the Jacobian of the transformation ( |3.36 ), ( 3.37 ) from the variable Q to 

{Qq, W} (the Jacobian does not contribute to the leading order correction calculated 
here, but is important for higher-order calculations |^5|, 0). Expanding up to the 
order W^, we get 

S-eff = 5o + - ]^{Sl) + + . . . . (3.40) 

The integral over the fast modes W can be calculated now using the Wick theorem 
with the contraction rules induced by the action Sw, 

{^iYW(Y)PW{Y')R) =n(r,r')(StrPStri?-StrPAStri?A) ; 
(Str[iy(r)P]Str[W^(r')P]) = n(r, rOStr(Pi? - PAPA), (3.41) 
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where P and R are arbitrary supermatrices. The diffusion propagator 11 is the solution 
of the diffusion equation 

- DV'U{r,, rs) = {nu)-'[S{r, - r^) - (3.42) 

with the Neumann boundary condition (normal derivative equal to zero at the sample 
boundary) and can be presented in the form 

n(r,r') = ^ E 7</>.(r)0.(r') , (3.43) 

where (^^{r) are the eigenfunctions of the diffusion operator —DV^ corresponding to 
the eigenvalues (equal to Dq^ for a rectangular geometry). As a result, we find 

= , 

(^2^ = 1 y d'^rdVn2(r, r')(5StrgoA + wStrgoA^)^ • (3.44) 
Substitution of Eq. j ^M) into Eq. (|^) yields 

2 

5'eff[Qo] = JsStrgoA + ^wStrgoAfc + A^(sStrgoA + wStrgoAA;)^ ; 
Zi Zi Lb 

A=^J d^rdVn^(r, r') = A E f-V = ^ ' (3-45) 

where is a numerical coefficient depending on the spatial dimensionality d and on 
the sample geometry. Assuming a cubic sample with hard-wall boundary conditions, 
we find 

1 °° 1 

ni, . . . ,nd = 
+ . . . + > 

yielding for = 1,2,3 the values ai = 1/90 ^ 0.0111, 02 ~ 0.0266, and 03 ~ 0.0527 
respectively.]^ In the case of a cubic sample with periodic boundary conditions we get 
instead 

1 °° 1 

= TTTTT 2^ \ 2^2 ' (3.47) 

(27r)4 ^ (nf + . . . + n^y 

ni,...,nd = -00 

nl + ... + nl> 

so that ai = 1/720 ~ 0.00139, 02 ~ 0.00387, and 03 ~ 0.0106. Note that for d < 4 the 
sum in Eqs. (|3.46| ), ( p.47| ) converges, so that no ultraviolet cutoff is needed. 

^^While speaking about d — 1, we mean a sample of a quasi-lD geometry, i.e. either a 2D strip 
b X L with ^ i or a 3D wire 61 x 62 x i with &i, 62 ^ L. For a strictly ID sample (a chain) the 
diffusive cr-model is not applicable. 
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Using now Eq.( p.35| ) and calculating the remaining integral over the supermatrix 
Qq, we finally get the following expression for the correlator to the order: 

^ sin^fvrs) . 9. ^ , 4arf 

R2is) = 1 + dis) - . \/ +Asin^ TTg ; A = ^. 3.48 



The last term in Eq. (|3.48|) represents the sought correction of order 1/ g"^ io the Wigner- 
Dyson distribution. 

The important feature of Eq.( p.48| ) is that it relates corrections to the smooth and 
oscillatory parts of the level correlation function, 

6R2{s) = A sin^ vrs = -(1 - cos 2tts) . (3.49) 
While appearing naturally in the framework of the supersymmetric cr-model, this fact is 



highly non-trivial from the point of view of semiclassical theory which represents 
the level structure factor K{t) (Fourier transform of R2{s)) in terms of a sum over 
periodic orbits. The smooth part of R2{s) corresponds then to the small-r behavior of 
K{t), which is related to the properties of short periodic orbits. On the other hand, 
the oscillatory part of R2{s) is related to the behavior of K{t) in the vicinity of the 
Heisenberg time r = 27r (t = 27r/A in dimensionful units), and thus to the properties 
of long periodic orbits. We will return to the non-universal corrections to K{t) below. 

The calculation presented above can be straightforwardly generalized to the other 
symmetry classes. The result can be presented in a form valid for all the three cases: 

Ri'\s) = (1 + A^.^) Ri%,Ms) , ^ = ^ , (3.50) 

where (3 = 1(2,4) for the orthogonal (unitary, symplectic) symmetry; -R2^rmt denotes 
the corresponding RMT correlation function (|2.71| )- (|2.73|) . 

For s ^ the RMT correlation functions have the following behavior: 

~ s ^ • 

-'^2, RMT — ^13^ 5 t> —>■ u , 

= C,^'-, C.= <^. (3.51) 

6 ' 3 ' 135 ^ ^ 

As is clear from Eq. (|3.50|) , the found correction does not change the exponent f3, but 
renormalizes the prefactor C/3: 

<M^(l+ (^ + ^)C'^+^ ^)c,,^ .^0 (3.52) 

The correction to C/j is positive, which means that the level repulsion becomes weaker. 
This is related to a tendency of eigenf unctions to localization with decreasing g. 

3.3.3 Stationary-point method. 

Let us return now to the behavior of the level correlation function in its high-frequency 
tail s g. We have already discussed the non-oscillatory part of R2{s) in this region. 
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see Eq. ( p.34|) . What is the fate of the oscillations in R2{s) in this regime? The 
answer to this question was given by Andreev and Altshuler who calculated i?2(s) 
using the stationary-point method for the cr-model integral ( |3.17| ) and treating the zero 
mode and the higher modes on equal footing. Their crucial observation was that on 
top of the trivial stationary point Q = A (expansion around which is just the usual 
perturbation theory), there exists another one, Q = kA, whose vicinity generates the 
oscillatory part of R2{s). (In the case of symplectic symmetry there exists an additional 
family of stationary points, see |^8[). The saddle-point approximation^ of Andreev and 



Altshuler is valid for s ^ 1; at 1 ^ s ^ (? it reproduces the above results of Ref. |^ 



(we remind that the method of [^] works for all s <^ g). The result of has the 
following form: 



(3.53) 
(3.54) 
(3.55) 



-^2,osc 




cos27rs ^. , 


-^2,080 




cos 2ns , , 


pSp 
-'i'2,osc 







where D{s) is the spectral determinant 

^w=^n(i+ 



COS Arts 
327r4 



2 A2\ -1 



D''(s] 



s^A 



(3.56) 



The product in Eq.( p.56D goes over the non-zero eigenvalues of the diffusion operator. 
This demonstrates again the relation between -R2,osc(s) and the perturbative part (|3.31|) , 
which can be also expressed through -D(s), 



Ri%is) - 1 



1 d'^ln D{s) 



(3.57) 



In the high-frequency region s 3> (7 the spectral determinant is found to have the 
following behavior: 



D(s] 



exp 



TT 



2sY'' 



(3.58) 



r(rf/2)c;sin(7rc//4) V ^ / 
so that the amplitude of the oscillations vanishes exponentially with s in this region. 



^^To avoid possible confusion, we remind that the matrices of the manifold (3.14) are exact saddle- 
points of the action (3.9) {i.e. they all have exactly the same action) for w = only. At finite cu 



this becomes a manifold of quasi-saddle-points, with the action difference determined by the term 
(— z7rs/2)StrQA, see Eqs. ( [2.61 ) and ( 3.16| ). The corresponding soft (cr-model) modes should be con- 
trasted to massive modes (describing fluctuations in Pi. 2, see Eq. ( 2.51 )), which have been integrated 
out in the course of derivation of the cr-model. Now, on this manifold of almost-saddle-points Q there 
are two true (even at non-zero lo) stationary points, namely, Q = A and Q — kA. In fact, we have 



already mentioned the existence of the second diagonal saddle point below Eq. (2.58). It is easy to 
see that the choice of signs S2 = — S4 = 1 produces there precisely the matrix Q = kA. 
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Taken together, the results of and ||4^ provide a complete description of the 



deviations of the level correlation function from universality in the metallic regime g 3> 
1. They show that in the whole region of frequencies these deviations are controlled by 
the classical (diffusion) operator governing the dynamics in the corresponding classical 
system. 

3.4 Spectral characteristics related to i?2(s). 
3.4.1 Spectral formfactor. 

The spectral formfactor is defined as the Fourier transform of the connected part 
R^\s) = R2{s) — 1 of the two-level correlation function, 

/oo 
4'^(s)e'""ds . (3.59) 
-oo 

By definition, K{t) (as well as R2{s)) is an even function, so that it is sufficient to 
discuss it at r > 0. In GUE it has the form 

Let us analyze, what kind of corrections to K2{t) is implied by the deviations of R2{s) 
from universality studied in Sec. For this purpose, let us use Eq. ( |3.49|) , forgetting 



for a moment about the condition of its validity {s g). The Fourier transformation 
of (ICTp then yields 



Att 

5K{t) = — [25(r) - 5(r - 2n) - 5(r + 27r)] . (3.61) 

Taking now into account the existence of the cutoff s ~ (7 for ( |3.49| ) leads to smearing of 
the (5-functions in (|3.61|) over an interval ~ 1/(7 around r = and r = ±27r, respectively. 
Thus, we conclude that 

5K{t) = 5Ko{t) + 5K2^{r) , r > , (3.62) 

where 5Kq{t) is located in the interval between r = and r ~ 1/(7, while 5K2-k{t) is 
concentrated in the interval of a width ~ 1/(7 around r = 27r. Since (5i?2(0) = 0, the 
integrals of 5Kq{t) and 5K2tt{j) are equal up to a sign, 

dr5ifo(T) = - j A.t6K2^{t) . (3.63) 

Furthermore, using the identity 

d^ 



we find 



dr5ir2.(r)~-j, (3.64) 
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so that the correction around the Heisenberg time r = 2-71 is negative. Since the decay 
of -R2,osc(s) at s > (7 is exponential, 5K2tt{t) is a smooth function. Using the fact that 
it is essentially located in an interval of width ~ 1/(7 and that A ~ l/(yf^, we conclude 
that the magnitude of 5K2t,{t) in the above interval is ~ \/g. The correction K2t,{t) 
has thus both the magnitude and the width of the order oilj g and leads to a rounding 
of the singularity in the spectral form-factor at the Heisenberg time r = 27r, as was first 



realized by Andreev and Altshuler |48|. In the quasi-lD case, the spectral determinant 



D{s), Eq. ( p.56| ), and thus the correction 5K2t,{t) can be calculated analytically, see 



As to 5Kq{t)^ its small-r behavior depends on spatial dimensionality, in view of the 
the Altshuler-Shklovskii "tail" ( |3.34|) . Taking the Fourier transform, we findQ 

5Ko{t) ~ -^{grf-'l' , r<l/g. (3.65) 

For semiclassical treatment of the spectral form-factor applicable for t/2'k ^ 1 [where 
it reproduces the formula ( ^.65| )] see Ref. PU| . 

3.4.2 Level number variance. 

Consider the variance ^2(3) of the number of energy levels in a spectral window of a 
width 5E = sA (the average number of levels in this interval being equal to s). It is 
easy to see that ^2(5) is related to the two-level correlation function as follows: 

1:2(3)= f dS{s -\s\)R^2\s) . (3.66) 

J — s 

The relation (|]6|) can be also presented in the following way 

^^1:2(8) = 1^ d~s R^^\~s) . (3.67) 

The level number variance is commonly used to characterize the long-range behavior of 
spectral correlations. The decay of the smooth partQ of the two- level correlation 
function at s 3> 1 in RMT implies the Ins behavior of S2(s). In particular, for GUE 
the large-s asymptotics reads 0, |^ 

S^rmt(s) = ^[1 + 7 + ln(27rs)] , s » 1 . (3.68) 



According to (|3.49|) , the correction to the RMT form of ^2(5) in a diffusive sample is 



small (and positive) at s <ti g and has the form 

Ay ( s 

-s^ = 2aA- 

2_ \g 

^^In 3D the spectral formfactor (|3.65| ) seems to diverge as r — > 0. It should be taken into account, 
however, that the above considerations are valid only for frequencies uj corresponding to diffusive mo- 



6E2{s) = = 2ad(-^ , s-^g . (3.69) 



tion, i.e. lu , where Tq is the elastic mean free path (denoted by t in Sec. 3.1). Correspondingly, 



the applicability of (3.65) is restricted by the condition r > t^A. 

^''The oscillatory part of i?2(s) is not important for the behavior of S2(s) at s ^ 1, because it gives 



a negligible contribution after integration ( 3.66 ) 
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On the other hand, at s ^ g the level number variance is determined by the Altshuler- 
Shklovskii behavior of R2{s) and is totally different from its RMT form: 



S ^ ' 



S2(s) ~ (^-J > S2,rmt(s) , s > (7 . (3.70) 

4 Eigenfunction statistics. 

Not only the energy levels statistics but also the statistical properties of wave functions 
are of considerable interest. In the case of nuclear spectra (the statistical description 
of which was the original motivation for the development of the RMT) , they determine 
fluctuations of widths and heights of the resonances |51|]. A more recent growth 
of interest to statistical properties of eigenf unctions in disordered and chaotic systems 
has been motivated, on the experimental side, by the possibility of fabrication of small 
systems (quantum dots) with well resolved electron energy levels p2| , |53| , 0, Fluc- 



tuations in the tunneling conductance of such a dot measured in recent experiments 



56, 57] are related to statistical properties of wavefunction amplitudes |58, 5£, 60, 31 
Furthermore, the eigenfunction fluctuations determine the statistics of matrix elements 
of the Coulomb interaction, which is important for understanding the properties of ex- 
citation and addition spectra of the dot It is also worth mentioning 
that the microwave cavity technique allows one to observe experimentally spatial fluc- 
tuations of the wave amplitude in chaotic and disordered cavities [|^, (though 
in this case one considers the intensity of a classical wave rather than of a quantum 
particle, all the results are equally applicable). 



4.1 Random matrix theory 

Within the RMT, the eigenfunction statistics has a very simple form in the limit ^ 1. 
The components of an eigenvector ip^^^ of a matrix H from the Gaussian ensemble 
become then uncorrelated Gaussian random numbers (real, complex, or quaternionic 
for P = 1, 2, and 4, respectively) with the distributions 

p(V^«)«exp|-^|V.«p} . (4.1) 

If we introduce the "intensity" 

y^:^=N\4^\' (4.2) 

(we have chosen the normalization (y) = 1), its distribution will then have the form 
of the x^-distribution with f3 degrees of freedom. In particular, for GUE and GOE we 
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pU(^) = , (4.3) 



e 



Equation ( [4.4|) is known as the Porter-Thomas distribution; it was originally introduced 
to describe fluctuations of widths and heights of resonances in nuclear spectra . 



4.2 Eigenfunction statistics in terms of the supersymmetric 
cr-model 

Theoretical study of the eigenfunction statistics in a dimensional disordered sys- 



tem is again possible with making use of the supersymmetry method ||69|, [7^, ^ . 
We consider the local intensity of an eigenfunction at some point fq of the sample, 
Z/i(i"o) = ^|V'r^(i"o)| (we have again normalized it to {y) = 1). In order to calculate the 
distribution function V{y), we introduce its moments, 

W = A{J2 \Uro)\''S{E - E,)) ^ m^o)?')E ■ (4.5) 

i 

The next step is to observe that /q(ro) is related for g = 2, 3, . . . to the following product 
of the Green's functions^ 

Ki,m{ro,ri) = {ro\{E-H + t7])-'\roy{ro\{E-H-zr])-'\ror 

^ [Gi(ro,ro)]'[Gf(ro,ro)r , (4.6) 

where /,m > 1 and l + m = q. Indeed, in the limit rj ^ (Ki^rn) possesses a singularity 
cx rf—^~^ determined by the probability oc to find an eigenvalue of H within a 
distance ~ 7] from E, in which case Ki^rn oc rj"'''"^. Extracting this singularity, we get 

^'('»' - ^'"•" 'V+''l'-2)f' ,1%(2'')'^"-'('^'..»('o. "» ■ (^■^) 

Note that at this stage we have for q > 2 the freedom in choosing I and m (constrained 
only by / + m = g). The final result will not, of course, depend on this choice. 

We express now the product of Green's functions as an integral over the supervector 
field 

KU^o,v) = TTT / D'^D'^^[Si{ro)Sl{ro)f[S2{ro)S;{ro)r 
llml J 

xexpj^ j d'^r'i>^r)L^/^{E + ir]A~ H)L^/^<^{r)~^ . (4.8) 



^^For GSE this consideration gives Vd/) = Aye~^'^. Note, however, that in terms of the electronic 
wave function this corresponds to defining y as the total (summed over the spin projections) intensity, 
y — + IV'iP)- For the distribution of the spin-projected intensity y — 2V\ipi\'^ the RMT result 



would have the same form (4.4) as for GUE. We do not consider the symplectic symmetry below; 
deviations from universality have in the symplectic symmetry class a form similar to the results for 
the systems of the unitary and the orthogonal symmetry. 
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For q — and q = I the moments are trivial, /o(ro) = 1, -^i(ro) = ^/V- 
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Proceeding in the same way as in the case of the level correlation function (Sec. |2.1|) , 
we represent the r.h.s. of Eq. ( [4.8| ) in terms of a a-model correlation function. In the 
case of the unitary symmetry the results reads 



^i,m(ro,^) 



-iTTu) 



m 
3 



^ '5iifbb('^o)Q^,bb(''o)Qi2,bb(''o)Q2i,bb(''o)e" 



S[Q] 



where S[Q] is the a- model action, 



S[Q] 



d'^r Str 



TTUD 



{VQY - nur]AQ 



(4.9) 



(4.10) 



We remind that the two pairs of indices of the Q-matrices refer to the retarded-advanced 
(1, 2) and the boson-fermion (b, f) decomposition, respectively. 

Since the preexponential factor in ( ^4.9| ) depends on the Q-field at the point fq only, 
it is convenient to introduce the function Y{Qq) as the result of integrating out all 
other degrees of freedom, 



Y{Qo) 



BQ{v)exp{-S[Q]} . 



(4.11) 



'Q{ro)=Qo 

With this definition, Eq. ( p.4|) takes the form of an integral over the single matrix Qq, 

I 



.l+m 



J 



m 
3 



^ Qo:ll,hhQo;2'2,hhQi;12,hhQi;21,hhYiQo) 



(4.12) 



For invariance reasons, the function Y{Qo) turns out to be dependent in the unitary 
symmetry case on the two eigenvalues 1 < Ai < oo and — 1 < A2 < 1 only, when 
the parametrization (|2.62|) , ( 2.63| ) for the matrix Qq is used. Moreover, in the limit 
— > (at a fixed value of the system volume, and thus of the level spacing A) only 
the dependence on Ai persists. 



(4.13) 



with relevant values of Ai being Ai ~ A/t] ^ 1. One can further show that in this 
asymptotic domain 

Qll,bbQ22,bb — Ql2,bb<521,bb , 

SO that all terms in in Eq. ( [4.12|) are equal. Evaluating the integral over all coor- 
dinates but Ai in the parametrization ( |2.62| ), (|2.63| ) of Qq, we get 



/,(ro) = ^g(g-l) J dzz^-'Y,{z). 



(4.14) 



Consequently, the distribution function of the eigenfunction intensity is given by |^ 
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(4.15) 



In the case of the orthogonal symmetry, Y{Qq) = Y{Xi, A2, A), where 1 < Ai, A2 < 
cxo and — 1 < A < 1. In the limit rj —>■ 0, the relevant region of values is Ai ^ A2, A, 
where 

Y{Qo)-^Y,{nlx^) . (4.16) 



The distribution of eigenfunction intensities is expressed in this case through the func- 
tion K as follows [EUi: 



1 rco i2 

2^2 d2 r / 



(4.17) 



In a metallic sample, typical configurations of the Q-field are nearly constant in 
space, so that the zero-mode approximation is expected to be a good starting point 
(see Sec. pl2| ). It amounts to approximating the functional integral ([4.9|), ( ^.11| ) by an 
integral over a single supermatrix Q, yielding 



YJz) 



(0,U) 



(4.18) 



which reproduces [after substitution in ( [4.15| ), (14.17]) ] the RMT results ( ^73|) and 
Therefore, like for the level correlations, the zero mode approximation yields the RMT 
results for the distribution of the eigenfunction amplitudes. To calculate deviations 
from RMT, one has to go beyond the zero-mode approximation and to evaluate the 
function Yi,{z) determined by Eqs. ( [4.11| ), ( |4.13| ) for a (i-dimensional diffusive system. 
In the case of a quasi- ID geometry this can be done exactly via the transfer-matrix 
method |7D[. The result depends crucially on the ratio L/C,, where L is the system 
size and ^ = 27ri'AD the localization length [A being the transverse cross-section of 
the sample). The two limiting cases correspond to the metallic regime L <^ C, (with 
the dimensionless conductance g = ^/L ^ 1) and to the strong localization regime 
L ^ C. For higher d, the exact solution is not available any more, and one should 
rely on approximate methods. Corrections to the "main body" of the distribution 
can be found by treating the non-zero modes perturbatively [J^, while the asymptotic 
"tail" can be found |]7T[ via an instanton method [j72|. For the quasi-lD geometry 
these approximate methods reproduced the results obtained earlier from the exact 
solution. 



4.3 Quasi-one-dimensional geometry 
4.3.1 Exact solution of the cr-model 

In the case of quasi-lD geometry an exact solution of the a-model is possible due to 
the transfer-matrix method. The idea of the method, quite general for one-dimensional 



problems, is in reducing the functional integral of the type ([4.11|) to solution of a differ- 
ential equation. This is completely analogous to constructing the Schrodinger equation 
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from the quantum-mechanical Feynman path integral. In the present case, the role of 
the time is played by the coordinate along the wire, while the role of the particle coor- 
dinate is played by the supermatrix Q. In general, at a finite value of the frequency rj in 
Eq. ( |4.1CI| ) (more precisely, 77 plays the role of imaginary frequency) , the corresponding 
differential equation is too complicated and cannot be solved analytically. However, a 
simplification appearing in the limit 77 0, when only the non-compact variable Ai 
survives, allows to find an analytical solution |3^, |6^, |70| of the ID cr-model0. 



There are several different microscopic models which can be mapped onto the ID 
supermatrix a-model. First of all, this is a model of a particle in a random potential 
(discussed above) in the case of a quasi-lD sample geometry. Then one can neglect the 
transverse variation of the Q-field in the cx-model action, thus reducing it to the ID 
form [|73], Secondly, the random banded matrix (REM) model (already mentioned 



in Sec. |2.5.4| ) has been mapped onto the ID cx-model ||3^, |6^, |7^. The REM model is 



relevant to various problems in the field of quantum chaos. In particular, the evolu- 
tion operator of a kicked rotor (paradigmatic model of a periodically driven quantum 
system) has a structure of a quasi-random banded matrix, which makes this system 



belong to the "quasi-lD universality class" (see Sec. for a more detailed discus 



sion). Finally, the lida-Weidenmiiller-Zuk random matrix model |]74| of the transport 



in a disordered wire can be also reduced to the ID cr-model, which allows to study 



analytically the wire conductance and its fluctuations ||7^, |75 |. 

The result for the function Ya.{y) determining the distribution of the eigenfunction 
intensity reads (for the unitary symmetry) 

YM = W^'\y^/L,T+)W^'\yaL,T_) . (4.19) 

Here ^ = 2711/0 A is the localization length (where A is the wire cross-section), and 
r_ = x/^, = {L — x)/^, with < x < L being the coordinate of the observation 
point ro along the sample. For the orthogonal sjTiimetry ^ is replaced by ^/2. The 
function W^^\z,t) satisfies the equation 

(4.0, 

and the boundary condition 

W^^\z,0) = l. (4.21) 

The solution to Eqs.( |4.20| ), ( ^.21| ) can be found in terms of the expansion in eigenfunc- 
tions of the operator z'^-^ — z. The functions 2z^^'^Ki^{2z^^'^), with K,y{x) being the 
modified Eessel function (Macdonald function), form the proper basis for such an ex- 
pansion [^], which is known as the Lebedev-Kontorovich expansion; the corresponding 
eigenvalues are — (1 + /^^)/4. The result is 

W^^\z,t)=2z'/^{k^{2z'/') + - r di2-^smh^KU2z'/^)e-'-^A . (4.22) 
I ^ Jo l + /i^ 2 J 

^^Let us stress that we consider a sample with the hard-waU (not periodic) boundary conditions in 
the longitudinal direction, i.e a wire with two ends (not a ring). 
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The formulas ( [4.15| ), (f4.17|) , ( [4.19|) , (f4.22| ) give therefore the exact solution for the 
eigenfunction statistics for arbitrary value of the parameter X = L/C,- The form of the 
distribution function is essentially different in the metallic regime X <^ 1 (in this 
case X = 1/g) and in the insulating one, X ^ 1. We discuss these two limiting cases 
below. 



4.3.2 Short wire 

In the case of a short wire, X = 1/g < 1, Eq. ( ^:T5|) , i ^A7\ ), j ^J^ , {^^^ yield 



p(u)(^) = [1 + ^(2 - 4y + y2) + . . .] 
exp{-y + fy2X + ...} ; 



= 7^ exp { 1 [-y + f y^X + . . .] } 



exp 



-2/3 



y < x-'i' 
y < x-'i' 
x~''' <y<x-' 
y 

y>x~'. 



(4.23) 

(4.24) 

(4.25) 
(4.26) 

(4.27) 



Here the coefficient a is equal to a = 2[1 — 3x(L — a;)/L^]. We see that there exist three 
different regimes of the behavior of the distribution function. For not too large ampli- 
tudes y, Eqs.( |4.23| ), ( [4. 241 ) are just the RMT results with relatively small corrections. 
In the intermediate range ( [4.25| ), ( [4.26| ) the correction in the exponent is small com- 
pared to the leading term but much larger than unity, so that V{y) ^ VRMxiy) though 
lnP(|/) ~ lnPijMr(?/)- Finally, in the large-amplitude region, (|4.271) , the distribution 
function V{y) differs completely from the RMT prediction.]^ 



4.3.3 Long wire 



In the limit of a long sample, X = L/{ > 1, Eqs. (P^ ), ( CT ), ( ^TTUQ , reduce 
to 



V 



u 



L 



(4.28) 
(4.29) 



with ^ = 2'KvAD as before, and u = |'?/'^(ro)|. Note that in this case it is not appropriate 
to use y = uV as a variable, since typical intensity of a localized wave function is 
u ~ 1/^^ in contrast to m ~ 1/V for a delocalized one. The asymptotic behavior of 
Eqs. ( [4.28| ), ( ^4.29] ) aX l/A^ has precisely the same form. 



V{u) ~ exp(-2/?v^^) 



(4.30) 



Note that Eq. ( 4.27 ) is not valid when the observation point is located close to the sample bound- 



ary, in which case the exponent of (4.27) becomes smaller by a factor of 2, see |77 
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as in the region of very large amplitude in the metallic sample, Eq. (|4.27|) . On this 
basis, it was conjectured in |7y] that the asymptotic behavior ( |4.27|) is controlled by the 
probability to have a quasi-localized eigenstate with an effective spatial extent much 
less than ^ ("anomalously localized state"). This conjecture was proven rigorously 
in |]7^ where the shape of the anomalously localized state (ALS) responsible for the 
large-M asymptotics was calculated via the transfer-matrix method. 

The transfer-matrix method allows to study, in the quasi-lD geometry, not only 
statistics of the eigenfunction amplitude in a given point, but also correlation functions 
of amplitudes in different points. Relegating a more extensive discussion of this issue 
to Sec. |4.5| , we mention here a distribution function which characterizes fluctuations 
in the rate of exponential decay of eigenfunctions (Lyapunov exponent). Specifically, 
let us consider the product of the eigenfunction intensity in the two points close to the 
opposite edges of the sample Xi — > 0, 0:2 ^ L, 



v = {2T,uDA^f\i;l{v,)i;l{v,)\ 



(4.31) 



The corresponding distribution function is found to be |T^, |T0 



V{-\nv) = F^^\-{l3\nv)/2X] 
,r^[(3-n)/2] 



/3 



1/2 



exp 



(3 (2X 



u- 



Tiu) 



8nX J " 1 8X V /3 

uT^[{l-u)/2] 



—— + Inf 



(4.32) 



Therefore, Inf is asymptotically distributed according to the Gaussian law with mean 
value (— Inw) = {2/f3)X = L/fi-KvAD and variance var(— Inf) = 2(— Inw). The same 
log-normal distribution is found for the conductance and for transmission coefficients 
of a quasi- ID sample from the Dorokhov-Mello-Pereyra-Kumar formalism |80[| . 

Note that the formula ( |4.32| ) is valid in the region o{v<^l (i.e. negative Inw) only, 
which contains almost all normalization of the distribution function. In the region 
of still higher values of v the log- normal form of 'P{v) changes into the much faster 
stretched-exponential fall-off oc exp{— 2\/2/3f ^/^}, as can be easily found from the exact 
solution given in The decay rate of all the moments (f*^), k > 1/2, is four 

times less than {—\nv) and does not depend on k: (f^) oc e~"^/^^. This is because 
the moments {v^), k > 1/2 are determined by the probability to find an "anomalously 
delocalized state" with t> ~ 1. 



4.4 Metallic regime (arbitrary d). 



For arbitrary dimensionality d, deviations from the RMT distribution V{y) for not too 
large y can be calculated via the method of Ref. described in Section p.3.2 
Applying this method to the moments ([4.7D, ( [4.9|) , one gets 



1+ 2^(9-1) + ••• 



(U) 



(4.33) 
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where 



K = n(ro,ro) = ^ 



True, 



Correspondingly, the correction to the distribution function reads 



l + -(2-% + y2) 



Similar results are obtained for the orthogonal symmetry class 

(2g-l)!! 



4(ro) 



■[l + nq{q-l) + .. 



(O) 



(4.34) 



(4.35) 



(4.36) 



1 + 



3y + 



y' 



+ 



(4.37) 



3 

2 V2 ' 2 

Numerical studies of the statistics of eigenfunction amplitudes in the weak-localization 
regime have been performed in Ref. [PT| for the 2D and in Ref. [P2| for the 3D case. 
The found deviations from RMT are well described by the above theoretical results. 
Experimentally, statistical properties of the eigenfunction intensity have been studied 



for microwaves in a disordered cavity For a weak disorder the found deviations 
are in good agreement with ( [4.37|) as well. 

As we see from the above formulas, the magnitude of the corrections is governed by 
the parameter k = n(ro, fq) (the one-diffuson loop in the diagrammatic language). In 
the quasi-one- dimensional case (with hard wall boundary conditions in the longitudinal 
direction), it is equal to 



K = n(ro,ro) 



9 



X 



1 



< X < L , 



(4.38) 



where x is the longitudinal coordinate of the observation point ro, so that Eqs. (|3.12|) , 
(|3.13|) agree with the results ( [4.23|) , ( [4.24| ) obtained from the exact solution. For the 
periodic boundary conditions in the longitudinal direction (a ring) we have k, = 1/Qg. 
In the case of 2D geometry, 

n(r,r) = — In^, (4.39) 
ng I 

with g = 27rz/D. Finally, in the 3D case the sum over the momenta n(r, r) = 
{ttuV)^^ X]q(-^^^)^^ diverges linearly at large q. The diffusion approximation is valid 
up to g ~ l^^; the corresponding cutoff gives n(r, r) ~ 1/2'kvDI = g^^{L/l). This di- 
vergency indicates that more accurate evaluation of n(r, r) requires taking into account 
also the contribution of the ballistic region {q > which depends on microscopic 
details of the random potential; see for details. 

The formulas ( 4.35| ), ( 4.37|) are valid in the region of not too large amplitudes, 
where the perturbative correction is smaller than the RMT term, i.e. at ?/ -C k"^/^. In 
the region of large amplitudes, y > k"^/^ the distribution function was found by Fal'ko 
and Efetov |7T| who applied to Eqs. ( |4.15| ), ( [4.17| ) the saddle-point method suggested 
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by Muzykantskii and Khmelnitskii [|72[. We relegate the discussion of the method to 
Sec. and only present the results here: 



z^"'/' < y < (4.40) 

(4.41) 



P(y) ~exp<' f + V + ...) } X 



r{y) ~exp<{ —^\n'^{Ky) 




Again, as in the quasi-one- dimensional case, there is an intermediate range where a 
correction in the exponent is large compared to unity, but small compared to the leading 
RMT term [Eq. ( [4.40| )] and a far asymptotic region ([4.41|) , where the decay of 'P{y) 
is much slower than in RMT. Similarly to the quasi-lD result ( |4.27|) , the asymptotic 
behavior (|4.41| ) is determined by anomalously localized states (see Q for a review). 

4.4.1 2D: Weak multifractality of eigenfunctions 

Since c? = 2 is the lower critical dimension for the Anderson localization problem, 
metallic 2D samples (with g 1) share many common properties with systems at 
the critical point of the metal-insulator transition. Although the localization length ^ 
in 2D is not infinite (as for truly critical systems), it is exponentially large, and the 
criticality takes place in the very broad range of the system size L <C ^. 

The criticality of eigenfunctions shows up via their multifractality. The multifractal 
structures first introduced by Mandelbrot |^ are characterized by an infinite set of 
critical exponents describing the scaling of the moments of a distribution of some 
quantity. Since then, this feature has been observed in various objects, such as the 
energy-dissipating set in turbulence [0, |8^, strange attractors in chaotic dynamical 
systems p^, B8|, p9|, BOl, and the growth probability distribution in diffusion-limited 



aggregation [91, 92, 93]; see Ref. |9^] for a review. 

The fact that an eigenfunction at the mobility edge has the multifractal structure 
was emphasized in |^5|] on the basis of renormalization group calculations done by 
Wegner several years earlier [PU]. For this problem, the probability distribution is 



just the eigenfunction intensity |'?/'^(r)| and the corresponding moments are the inverse 
participation ratios (IPR's), 



Pn 



2q, 



(4.42) 



The multifractality is characterized by the anomalous scaling of Pq with the system 
size L, 

Pq oc L-^'^"-^) = (4.43) 

with Dq different from the spatial dimensionality d and dependent on q. Equivalently, 
the eigenfunctions are characterized by the singularity spectrum f{a) describing the 
measure L^'^"^ of the set of those points r where the eigenfunction takes the value 
|'?/'^(r)| oc L~°'. The two sets of exponents r(g) and /(a) are related via the Legendre 
transformation, 

r(g) = qa- f{a) ; f{a) = q ; r'(g) = a . (4.44) 
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100, 101 



By now, the multifractality of critical wave functions is confirmed by numerical simu- 



lations |9^, |9^, |99[ ; for more recent reviews see Refs 

We turn now to the wave function statistics in 2D. We note first that the formulas 
(|]33|), (|]36|) for the IPR's with q < can be rewritten in the 2D case (with (|]3|) 

taken into account) as 

pRMT 

where is the RMT value of equal to glL-^fa-i) for GUE and (2g - l)!!L-2('/-i) 

for GOE. We see that ( 4.45 ) has precisely the form (^4.43 ) with 



(4.45) 



2 - 



(4.46) 



As was found in [|7T|, the eigenfunction amphtude distribution ( |4.40| ), ( [4.41| ) leads to 
the same result ( ^4.46|) for all q <C '^/Sng. Since the deviation of Dq from the normal 
dimension 2 is proportional to the small parameter l/irg, it can be termed "weak 
multifractality" (in analogy with weak localization). The result ( |4.46| ) was in fact 
obtained for the first time by Wegner |]5B| via the renormalization group calculations. 

The limits of validity of Eq. (^4.46|) are not unambiguous and should be commented 
here. The singularity spectrum f{a) corresponding to ( |4.46|) has the form 



/(«) 



pTcg 



a 



so that f{a± = 0) for 



1 ± 



1 



(4.47) 



(4.48) 



{2f37rgy/\ 

If a lies outside the interval (a_, a+), the corresponding /(a) < 0, which means that 
the most likely the singularity a will not be found for a given eigenfunction. However, 
if one considers the average (Pg) over a sufficiently large ensemble of eigenfunctions, 
a negative value of /(a) makes sense (see a related discussion in ||103| |). This is the 
definition which was assumed in |7l| where Eq. (|4.46| ) was obtained for all positive 
q < 2(37rg. 

In contrast, if one studies a typical value of Pg, the regions a > a+ and a < a_ will 
not contribute. In this case, Eq. ( |4.46| ) is valid only within the interval q^ < q < q+ 
with q± = ±{2f3Tcgy^'^; outside this region one finds ||104| , |105 



T{q)^Dg{q-l] 



qa^ 
qa+ 



q> q^ 
q<q^ 



(4.49) 



Therefore, within this definition the multifractal dimensions Dg saturate at the values 
a+ and «_ for q — >■ +oo and q — ^ — oo respectively. This is in agreement with results 
of numerical simulations P7| , |100| , |101|| . 
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4.5 Spatial correlations of eigenfunction amplitudes. 

Correlations of amplitudes of an eigenfunction in different spatial points are character- 
ized by a set of correlation functions (we consider, as usual, the unitary symmetry for 
definiteness) 

A{J2r^irl)Ur[)...r^ir,)Ur',nE-E,)) ^ (^*(ri)^(r;) . . . ^*(r,)^(r;)) . (4.50) 

i 

Using the supersymmetry formalism and performing the same transformations that 
have led us to Eq. (|]9|), we get 

= - 2V(,'-l)! i^i"o(^""^)"' / DQEri;^^-M..bb(r.,r;(,)e-^[«] , (4.51) 

^ ' a i=l 

where the summation goes over all transpositions a of the set {1, 2, . . . , g}, pj is equal 
to 1 for i = 1, . . . , g — 1 and to 2 for i = q, and g is the Green's function in the field 
Q(r), 



E- — + i—] . (4.52) 

Taking into account that the field Q{r) varies only weakly on the scale of the mean 
free path / yields^ 

g{r,, ra) ~ ReGA(ri - r^) - t ImGA(ri - r2)g(ri) , (4.53) 

For |ri — ^ / the Green's function (7(ri,r2) vanishes exponentially, in view of 
G'a(i") oc e^**/^'. Since in the limit 77 ^ the characteristic magnitude of Q is oc 
A/?7 S> 1 [see the text around Eq. (f4.13|) ], the real part ReG'A(i") can be neglected in 
all the Green's function factors in ([4.51|) , and only the products of the imaginary parts 
survive. The imaginary part ImGA(r) is given explicitly by 



I"^^A(r) _ r -r/2i ^ j Mpft) , 2D 



ppr 

Let us note that /f('") is determined by microscopic (short-scale) dimensionality of the 
sample rather than by its global geometry. In particular, a quasi-lD sample may be 
microscopically of 2D (strip) or 3D (wire) nature. 

^^Since Q is a slowly varying field, the argument of Q in the second term of ( 4.53| ) can be chosen 
to be either ri or r2, or (ri + r2)/2. 
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4.5.1 Zero-mode approximation. 

In the zero-mode approximation, Eq. ( [4. 511 ) reduces to 



(T 1 = 1 

Equation ( [4.56 ) implies that the wave function ip{r) has a global Gaussian distribution, 



P{^(r)} ocexpj-^ j d'^rdVV^*(r)ir(r,r')^(r') 



(4.57) 



determined by the correlation function 

V{r{r)ij{r')) = M\r-r'\) , (4.58) 

the kernel K{r, r') being the operator inverse of V^^f^dr — r'\). An analogous consid- 
eration for the orthogonal symmetry class (when ip{r) is real) leads, in the zero- mode 
approximation, to the same conclusion. The result ([4.57|) , ([4.58|) was first obtained by 



Berry | 106 | from the conjecture that a wave function in a classically chaotic system is 
given by a random superposition of plane waves.0 

As has been explained above, quite generally there are corrections to the zero- 
mode approximation induced by the diffusion modes. They change the eigenfunction 
correlations qualitatively by inducing correlations on long scales r ^ / (which are 
exponentially small in the zero-mode approximation). Similarly to our discussion of 
the wave function fluctuations (Sec. ^^), we proceed by first presenting results of 
the exact solution in the quasi-lD case and then consider the metallic regime for an 
arbitrary dimensionality. 

4.5.2 Quasi-ID geometry. 



In the case of the quasi-lD geometry of the sample the method described in Sec. ^^3 
allows to calculate analytically all the multipoint correlation functions ( [4.50| ); the re 
suits are presented in terms of multiple integrals of the type ( [4.22|) , see reviews ^0 



Without going into technical details, we quote here some important conclusions con- 
cerning the global statistics of eigenf unctions. A wave function ipi^) can be represented 
as a product 

ip{r) = $(r)^(x) , (4.59) 

where x is the coordinate of the point r along the sample. Here $(r) is a quickly 
fluctuating in space function, which has the Gaussian statistics ( [4.57| ), ( [4.58|) obtained 



30 



Historically, the equivalence of the Berry's conjecture and the zero-mode supersymmetry calcu- 



lation has been established in a somewhat convoluted and less general way. Prigodin et al |107| 
calculated, in the zero-mode approximation for the cr-modcl, the joint distribution function V{u,v) 
of the wave function amplitudes in two different spatial points, u — {^"^{r)], v = |-0^(r')|. Srednicki 



then demonstrated |10S| that the same results for V{u,v) are obtained from the Berry's conjecture of 
Gaussian statistics of ipir). We have demonstrated this equivalence above in a more transparent and 
general form. 



45 



above from the OD cr-model. The (real) function \I/(x) determines, in contrast, a 
smooth envelope of the wave function. Its fluctuations are long-range correlated and 
are described by the probability density 



P{^(a;)}oce-^/2/5«e-[^(°)+^(^)]/2gxp|-^ej^ dx | ^ (^"' Jdxe'-l^ , 

(4.60) 

where 6{x) = ln\&^(a;). 

The physics of these results is as follows. The short-range fluctuations of the wave 
function (described by the function $(r)) originate, as explained above, from the su- 
perposition of plane waves with random amplitudes and phases leading to the Gaussian 
fluctuations of eigenfunctions with the correlation function ( [4.55|) . The second factor 
\l/(x) in the decomposition ( [4.59|) describes the smooth envelope of the eigenfunction 



(changing on a scale ^ /), whose statistics given by ([4.60|) is determined by diffusion 
and localization effects. This factor is responsible for the long-range correlations {i.e. 
those on scales ^ I) of the wave amplitude. In the metallic regime, C,/L = g ^ 1, 
\l/(x) fluctuates relatively weakly around unity and the above long-range correlations 
of ijj{r) are parametrically small (see Sec. |4.5.3| ). In the opposite regime of a long wire. 



L/^ S> 1, the strong localization manifests itself in extremely strong spatial correlations 
of eigenfunctions. 

Finally, we compare the eigenfunction statistics in the quasi- ID case with that in a 
strictly ID disordered system. In the latter case, the eigenfunction can be written as 

ipm{x) = )J^cos{kx + (5)^(a;) , (4.61) 

where \l/(a;) is again a smooth envelope function. The local statistics of ipiBix) (i.e. 
the moments ( |4.5| )) was studied in | p.09|| , while the global statistics (the correlation 
functions of the type ( 4.501) ) in | 110 |. Comparison of the results for the quasi-lD 



and ID systems shows that the statistics of the smooth envelopes \l/ is exactly the 
same in the two cases, for a given value of the ratio of the system length L to the 
localization length (equal to jSiruAD in quasi-lD and to the mean free path / in ID). 
The equivalence of the statistics of the eigenfunction envelopes implies, in particular, 
that the distribution of the inverse participation ratios Pq [see Eq. ( |4.42| )] is identical in 
the ID and quasi-lD cases. This is confirmed by explicit calculations of the distribution 
function of P2, see Refs. 1TT1|, |TT2[ (ID) and Refs. [|TT3l ^ (quasi-lD). 



4.5.3 Metallic regime (arbitrary d). 

Correlations of eigenfunction amplitudes in the regime of a good conductor can be 



again studied via the method of Ref. described in Sec. |3.3.2| ; see Refs. |114 
The result has the form of the expansion in powers of the diffusion propagator 11. In 
particular, for the simplest correlation function showing long-range correlations we find 
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(up to the linear-in-II terms) 



V^{|*(r0ri*(r2)n 



1-2 



l + -n(ri,rO 



+ ^n(ri,r2) 



(4.62) 



The last term on the r.h.s. of ([4.62|) describes the long-range correlations between 
|'?/'(ri)p and |'?/'(r2)P induced by diffusion (or, in other words, by classical dynamics). 

In a similar way one can calculate also higher order correlation functions of the 
eigenfunction amplitudes. In particular, the correlation function (|'?/'''(ri)||'?/'^(r2)|) de- 
termines fluctuations of the inverse participation ratio P2, the result for the relative 
variance of 5(-P2) = var(P2)/ (-^2)^ being 



m) = I / 



d^'rd'^r' 



32arf 



(4.63) 



with the numerical coefficient ad defined in Sec. |3.3.2| (see Eqs. (|3.46| ), ( p.47|) ). The 
fluctuations ([4.63|) have the same relative magnitude (~ 1/(7) as the famous universal 
conductance fluctuations. Note also that extrapolating Eq. ([4.63|) to the Anderson 
transition point, where (7 ~ 1, we find 6{P2) ~ 1, so that the magnitude of IPR 
fluctuations is of the order of its mean value [which is, in turn, much larger than in the 
metallic regime; see Eq. ( |4.43| )]. 

Equation (WM) can be generalized onto higher IPR's Pq with q > 2, 



var(P,) 2 



rg'(g-i)' 



d'^rd'^r' 



r, r 



8g2(g-l)2a^ 



(4.64) 



so that the relative magnitude of the fluctuations of Pq is ~ g(g 
the higher irreducible moments (cumulants) ((P")), n = 2, 3, . . ., 



- V)/g. Furthermore, 
have the form 



[n 




1)! 




2 




{n 




1)! 




2 




{n 




1)! 


2 



d^Vi . . . d^Vn 



n(ri,r2)...n(r„,ri) 



Tr 



P 



g(g-l)n 



1) 



i- 



(4.65) 



where 11 is the integral operator with the kernel n(r,r')/V". This is valid provided 
q'^n <^ 2j3'n-g. Prigodin and Altshuler | 115 | obtained Eq. (|4.65 ) starting from the 
assumption that the eigenfunction statistics is described by the Liouville theory. Ac- 
cording to ( [4.65|) , the "central body" of the distribution function V{Pq) of the IPR Pq 
(with q'^/Png <^ 1) is determined by the spectrum of eigenvalues of the diffusion 
operator —DV^. For a more detailed study of this distribution function see [[11 5| , ^. 

The perturbative calculations show that the cumulants of the IPR's are correctly 
reproduced (in the leading order in 1/g) if one assumes [|115|| that the statistics of the 
eigenfunction envelopes \ip'^ {r)\smooth = e^*^'"^ is governed by the Liouville theory (see 



e.g. ||116| , |117|| and references therein) defined by the functional integral 



—e - 1 ) exp 



(4.66) 
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The "tails" of the IPR distribution function governed by rare reahzations of disorder 
are described by saddle-point solutions which can also be obtained from the Liouville 
theory description ( |4.66| ), see 0. The multifractal dimensions ( [4.46| ) in 2D can be 



reproduced by starting from the Liouville theory as well ||116| , |117|| . It should be 



stressed, however, that this agreement between the supermatrix cr-model governing 
the eigenf unctions statistics and the Liouville theory is not exact, but only holds in 
the leading order in 1/ g. In particular, the Liouville theory does not describe the wave 
function localization by weak disorder in 2D. 

Up to now, we have considered correlations between amplitudes of one and the same 
eigenfunction at different spatial points. One can also study correlations of different 



eigenfunctions, see [114]. Understanding of both types of correlations is important 
for evaluation of fluctuations of matrix elements {e.g. those of Coulomb interaction) 
computed on eigenfunctions ipk of the one-particle Hamiltonian in a random potential. 
Such a problem naturally arises when one investigates the effect of interaction onto 



statistical properties of excitations in a mesoscopic sample (see Refs. [|^, |6J, |63[|). 



4.6 Anomalously localized states and long-time relaxation. 

In this subsection we discuss one more method that can be used within the superma- 
trix cr-model formalism to investigate statistical properties of eigenfunctions. This is 



the instanton method introduced by Muzykantskii and Khmelnitskii in order to 
calculate the long-time dispersion of the average conductance G{t). Soon after the pa- 
per [|72[ appeared, it was realized that the method allows one to study the asymptotic 
behavior of distribution functions of different quantities, including relaxation times, 
eigenfunction intensities, local density of states, inverse participation ratio, and level 
curvatures. These asymptotics are determined by rare realization of disorder producing 
the states which show much stronger localization features than typical states in the 



system - anomalously localized states already mentioned in Sec. [4.3.3| and |4.4| . For a 



review of the obtained results and relevant references the reader is referred to fT^, §]• 
In fact, we have already quoted the results obtained in this way by Fal'ko and Efetov 
7l| for the "tail" of the eigenfunction statistics, see Eqs. ( [4.40|) , ([4.41|) . Here we wish 



to present the ideas of the method by discussing the original problem considered by 
Muzykantskii and Khmelnitskii. 

Let us consider (following Refs. |72|, |118|| ) the asymptotic (long-time) behavior of 



the relaxation processes in an open disordered conductor. One possible formulation 
of the problem is to consider the time- dependence of the average conductance G{t) 
defined by the non-local (in time) current-voltage relation 

lit) = I dt'G{t - t')V{t') (4.67) 

J — oo 

Alternatively, one can study the decay law, i.e. the survival probability Ps{t) for a 
particle injected into the sample at t = to be found there after a time t. Classically, 
Ps{t) decays according to the exponential law, Ps{t) ~ e"*/*°, where is the lowest 
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eigenvalue of the diffuson operator — DV^ with the proper boundary conditions. The 
time has the meaning of the time of diffusion through the sample, and t^^ is of the 
order of the Thouless energy. The same exponential decay holds for the conductance 
G{t), where it is induced by the weak-locahzation correction. The quantities of interest 
can be expressed in the form of the a-model correlation function 

G{t), m ~ I ^e--*y" DQ(r)A{Q}e^^[«] , (4.68) 

where S[Q] is given by Eq. ( p.l6| ). The preexponential factor depends on the 

specific formulation of the problem but is not important for the leading exponential 
behavior studied here. 

The main idea of the method is that the asymptotic behavior is determined by a 
non-homogeneous {i.e. r-dependent) stationary point of the action S[Q] (instanton). 
Such a stationary point is found by varying the exponent in Eq. (|4.68|) with respect to 
Q and ur, which yields the equations 



2DV(gvg) + iu;[A, Q] = (4.69) 



T 



j d'^rStr(Ag) =t (4.70) 



[We assume unitary symmetry; in the orthogonal symmetry case the calculation is 
applicable with minor modifications.] It remains 

i) to find a solution of Eq.( [4.69|) (which will depend on u); 

ii) to substitute it into the self-consistency equation ([4.7CI|) and thus to fix as a 
function of t; 

iii) to substitute the found solution Qt into Eq. (|4.68| ), 

m ~ exp I^Str I (VQ,)^! . (4.71) 

Note that Eq. ( ^.69| ) is to be supplemented by the boundary conditions Q = A at 
the open part of the boundary {i.e. boundary with an ideal metal) and (|3.28| ) at the 
insulating part of the boundary (if it exists). 

It is not difficult to show that the solution of Eq. (|4.69|) has in the standard 
parametrization ( p.62|) , (|2.63|) the only non-trivial variable - bosonic eigenvalue Ai = 
cosh 6*1, all other coordinates being equal to zero. As a result, Eq. ([4.69|) reduces to an 
equation for ^i(r) (we drop the subscript "1" below): 

V^^ + ^sinh^ = , (4.72) 



the self-consistency condition ( [4.70|) takes the form 

■ j d'^r(cosh^- 1) =t , (4.73) 



TTZ/ 
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and Eq. ([4 .711) can be rewritten as 

lnP,(t) = -^ I d'riVey . (4.74) 

For sufficiently small times, 9 is small according to ( |4.73| ), so that Eqs. ([4 .721) , ( |4.73| ) 
can be linearized: 

+ 27^ = ; 27 = iu/D ; (4.75) 



TTZ/ 



This yields 



d\e' = t . (4.76) 



2t \ 

— Mr), (4.77) 



where 0i is the eigenfunction of the diffusion operator corresponding to the lowest 
eigenvalue ei = t^^. The survival probability ( [4.74|) reduces thus to 



In m = ^ / eV'e = -^ld're' = -t/t^ , (4.78) 



as expected. Eq. ( |4.78|) is valid (up to relatively small corrections) as long as 6* ^ 1, i.e. 



for tA ^ 1. To find the behavior at t > as well as the corrections at t < A^^, one 
should consider the exact (non-linear) equation ( [4.72| ), the solution of which depends 
on the sample geometry. 

4.6.1 Quasi-ID geometry 

We consider a wire of length L and a cross-section A with open boundary conditions 
at both edges, 0{-L/2) = 0{L/2) = 0. Equations ( ^T2|) , ( ^73|) take the form 



rL/2 

e" + 2-fsmhe = , / dx(cosh^ - 1) = t/vrz/A . (4.79) 

J-L/2 



r.L/2 

The solution of ( [4.79|) yields the log- normal asymptotic behavior of Ps{t) at large times 



lnP,(t) ~ -^ln2(tA) ; tA > 1 , (4.80) 

with g = 2711' AD /L ^ 1 being the dimensionless conductance. 

Equation ( [4. 801 ) has essentially the same form as the asymptotic formula for G{t) 
found by Altshuler and Prigodin ||119|| for a strictly ID sample with a length much 



exceeding the localization length: 



G{t) r^exp^--\n\t/r)j (4.81) 
If we replace in Eq. ( [1.81| ) the ID localization length / by the quasi-lD localization 



length PnuAD, we reproduce the asymptotics ( [4.801) . This is one more manifestation 
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of the equivalence of statistical properties of smooth envelopes of the wave functions in 
ID and quasi-lD samples (see Sec. f4.5.2|) . Furthermore, the agreement of the results 
for the metallic and the insulating samples demonstrates clearly that the asymptotic 
"tail" ( [4. 801 ) in a metallic sample is indeed due to anomalously localized eigenstates. 



As another manifestation of this fact, Eq. ( |4.80| ) can be represented as a superposi- 
tion of the simple relaxation processes with mesoscopically distributed relaxation times 

Ps{t) ~ J dt^e-'/'^Vm) (4.82) 
The distribution function V{t^) then behaves as follows: 

P(t<^) ~exp|-^ln2((7At^)| ; > -L = . (4.83) 

Since the Thouless energy tj^^ determines the typical width of a level of an open system, 
the formula ( [4.83|) concerns indeed the states with anomalously small energy widths 



The saddle-point method allows us also to find the corrections to Eq. ( [4.78| ) in the 
intermediate region tD<t< |T2T|, i. 



lnP,{t) = ^(l--^^ + ...] , (4.84) 



with ti) = L'^/tc'^D. Equation ( |4.84| ) is completely analogous to the formula ( |4.25|) 



( |4.26|) for the statistics of eigenfunction amplitudes. It shows that the correction to 
the leading term — t/to in ln-Pg becomes large compared to unity at t > tD^/g, though 
it remains small compared to the leading term up to t ~ gtD ~ A~^. The result ( [4.84| ) 



was also obtained by Frahm ||122|| from rather involved calculations based on the equiv- 



alence between the ID cr-model and the Fokker-Planck approach and employing the 
approximate solution of the Dorokhov-Mello-Pereyra-Kumar equations in the metallic 
limit. The fact that the logarithm of the quantum decay probability, InPs(t), starts 
to deviate strongly (compared to unity) from the classical law. In Pf{t) = —t/t-Q at 
i ^ ^D\/^ was observed in numerical simulations by Casati, Maspero, and Shepelyansky 
123|| . For related results in the framework of a random matrix model see Sec. |4.6.3| . 



4.6.2 2D geometry 

Considering a 2D disk-shaped sample of a radius i?, one gets the following long-time 
asymptotics of P^{t) (or G{t)) l7|, |lT8[: 



P,{t) ~ {t/\)-^^^s , 1< tA < {Rllf (4.85) 

Ps(t) - exp {-^^^} , tA » {R/lf (4.86) 

where g = 2'kuD is the dimensionless conductance per square in 2D and r is the mean 
free time. Equivalently, these results can be represented in terms of the distribution 
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function ^(t^) of relaxation times, 



- { i: r .,,^nHU/rU . (4.87) 




4 ln(iJ/Z) 



where — R^/D is the time of diffusion through the sample. 

The far log-normal asymptotics ( [4.86| ) agrees with the one obtained earlier by Alt- 



shuler, Kravtsov, and Lerner from the renormalization-group (RG) treatment. Anal- 
ogous agreement between the instanton and the RG calculations was found for the 



asymptotics of the local DOS distribution function ||124|| . Note that the instanton 
method is superior to the RG treatment in several respects: (i) it is not restricted to 
2D or 2 + e dimensions; (ii) it is much more transparent physically, since the stationary- 
point solution 6{r) describes directly the spatial shape of the anomalously localized 
state, iV'^(r) Ismooth oc e^^^^; (iii) in some cases it allows to find intermediate asymptotics 
[see e.g. Eq. ( [4. 851 )] missed by the RG calculation. 



4.6.3 Random matrix model 

Here we mention briefly the results on the quantum decay law obtained by Savin and 



Sokolov ||125|| within the RMT model. This will allow us to see the similarities and the 
differences between the diffusive systems and the random matrix model. The model 
describes a Hamiltonian of an open chaotic system by a Gaussian random matrix 
coupled to M external (decay) channels. The found decay law has the form 

P,{t) ~ (1 + Tt/M)-^ , (4.88) 

where F = MT A/27r is a typical width of the eigenstate, with T characterizing the 
channel coupling (T = 1 for ideal coupling). In this case, the product MT plays the 
role of the dimensionless conductance g (in contrast to the diffusive case where g is 
governed by the bulk of the system, here it is determined by the number of decay 
channels and the strength of their coupling). For not too large t {tAT ^1), Eq. ( [4.88| ) 
yields the classical decay law, Ps{t) ~ e~*'", with corrections of the form 

InPs(t) = -tF(l-Ft/2M + ...) , (4.89) 



which is similar to the results found for the diffusive systems (see Eq. (|4.84| )). At large 
t ^ (AT)~^, the decay takes the power-law asymptotic form ||126| 



InPs(t) ~ -Mln(Ft/M) , (4.90) 
which is to be compared with Eqs. (|4.80| ) and ([4.85|) , (|4.86|) . 



5 Supersymmetry approach to the Quantum Chaos. 

The aim of this section is to review recently emerged ideas concerning application of 
the supersymmetry formalism to chaotic ballistic systems ("billiards"). It should be 
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stressed that this field is quite young, and most of the calculations performed so far 
are less rigorous than in the diffusive case; the limits of applicability of the obtained 
results are not well understood yet. In this sense the status of this section is somewhat 
different from the preceding part of this lecture course. 



5.1 Introduction: What have we learned from the diffusive 
problem. 

It is instructive to begin by summarizing what we have learned concerning the diffusive 
problem. As has been explained in Sec. ^, the statistical properties of energy levels 
and eigenfunctions in a diffusive disordered sample are described by the diffusive su- 
permatrix cr- model (|3.16|) . In the metallic regime, when the dimensionless conductance 
is large, g ^ 1, and the eigenfunctions are not localized {i.e. cover roughly uniformly 
the whole sample volume), the zero- mode approximation is a good starting point for 
treating the cx-model. It reduces the problem to the OD cr-model, yielding the RMT 
results for the level and eigenfunction statistics. Deviations from RMT (which are 
typically small in the "body" of the distribution functions but become large in the 
"tail" ) are controlled by the diffusion modes. In contrast to the universal RMT results, 
these deviations are system-specific, since they depend on the sample dimensionality, 
shape and size, as well as on the disorder strength. More specifically, the deviations 
are controlled by the diffusion operator —DV^, which determines the quadratic part 
of the action of the diffusion modes, 

S[W] = Y / d"^!- Striy2i[-^V2 - iuj]Wi2 . (5.1) 

The magnitude of the deviations from RMT is controlled by the small parameter g'^ ^ 
A/ei, where ei is the lowest non-zero eigenvalue of the diffusion operator (the Thouless 
energy) . 



5.2 Ballistic cr-model. 

Let us now turn to the case of ballistic chaotic systems, i.e. to the quantum chaos. 
We will first guess what the field-theoretical description of such a problem should look 
like and then will discuss how it can be derived. 



5.2.1 Heuristic arguments. 

A bulk of numerical simulations data have unambiguously demonstrated that gener- 
ically the statistics of energy levels and eigenfunctions amplitudes in a ballistic sys- 
tem (billiard) whose classical counterpart is chaotic is well described by RMT (this is 



known as Bohigas-Giannoni-Schmit conjecture [127]) . Since we already know that the 



OD (T-model is essentially equivalent to RMT, we expect that the sought field theory 
of quantum chaos should reduce in the leading (zero-mode) approximation to the OD 



53 



cr-model described in Sec. p.4[ Therefore, the field variable of this theory should be the 
supermatrix field belonging to the same coset space as the field Q{r) of the a-model. 

Furthermore, in analogy with the diffusive case, we expect deviations from RMT 
to be controlled by the modes of density relaxation in the system. Indeed, deviations 
from ergodicity in a diffusive sample are physically due to the fact that the process of 
a particle spreading over the whole sample volume is not instantaneous but requires 
a time ~ t-^, the relaxation being governed by the diffusion operator. Since the mo- 
mentum relaxation in a diffusive sample takes place on a much shorter time scale (of 
the order of the mean free time r), the diffusion operator (describing the dynamics on 
time scales ^ r) is an operator in r-space only. This is why the Q-field depends on 
the coordinate but not on the velocity in the diffusive case. Such a separation of the 
fast and slow dynamics is not applicable to a generic clean sample. Therefore, in this 
case the supermatrix field Q should depend both on the coordinate r and the velocity 
direction n = v/vp, Q = Q{r, n) (the absolute value of the velocity |v| = vp being fixed 
by the energy conservation). The classical dynamics in the phase space is governed by 
the Liouville operator 

C = {.,H} = —--——, (5.2) 
dp dr dr dp 

which reduces for the case of a billiard (no potential energy inside) to 

d 

C = fpn— (5.3) 

OT 

supplemented by the boundary conditions corresponding to the particle refiection by 
the sample boundary. It is clear from what has been said above that the Liouville 
operator £ is expected to replace the diffusion operator in the cr-model description. 
Therefore, the analog of ( |5.1|) should read^ 



(5.4) 



S[W{v,n)] = y/ d'^rdnStriy2i(r,n)[£-zc<;]W^i2(r,n) . 

It remains to restore the action in terms of the Q-field from the quadratic form 
It turns out, however, that in the ballistic case the corresponding action cannot be 
written in a simple form in terms of the Q-field,^ and it is more convenient to write 
it in terms of the T-matrix field parametrizing the cx-model manifold according to [cf. 

(PI), (iH)] 

g(r, n) = T(r, n)Ar"^(r, n) . (5.5) 

Using that T = 1 — W/2 + . . ., that C is the first-order differential operator and 
that the action should be invariant at — with respect to the global rotations 
Q{r,n) UQ{r,n)U~^ {i.e. with respect to T(r,n) — > [/T(r,n)), one concludes that 
the only allowed form is 

S[Q] = Y / ^"^^^^ Str[-2T~i(r, n)CT{r, n)A - tuQ{r, n)A] . (5.6) 

^^The angular integral J dn ... is assumed to be normalized to unity, J dn = 1. 
^^The action can be expressed in terms of the Q-ficld only at the expense of introduction of an 
additional coordinate, with the action taking the Wess-Zumino-Witten form, see [128|. 
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5.2.2 Ballistic cr-model from disorder averaging. 

The ballistic cr-model action was derived for the first time by Muzykantskii and Khmel- 
nitskii (MK) | ]1 28| | . Their starting point was a disordered system, and they followed 



the route outlined in Sec. ^71] up to Eq. (|3.9|) . Their further aim was to derive a theory 



which describes the physics at all momenta q <^ kp and not only at g ^ i.e. which 
is quasiclassical but is not restricted to the diffusion approximation. Clearly, in this 
case one cannot use the gradient expansion leading to the diffusive cr-model (|3.16|) . 
Instead of this, MK employed an analogy with the Eilenberger quasiclassical approach 
in the kinetic theory of disordered superconductors. Along these lines, they defined the 
field Q{r,n) as a result of application of the following two operations to the Green's 
function g{r,r') [defined in Eq. ( |3.1(]| )]: 

(i) the Wigner transformation 

(ii) integration over the kinetic energy ^ = frdpl —Pf), 

Q(r,n) = - / d^g{r,n-^) . (5.8) 
n J vf 

They were then able to derive the ballistic cr-modelQ ( |5.6| ) with an additional term 
describing the scattering by impurities, 

^imp[Q] = X / d''rdncin't/7(r; n, n')Strg(r, n)Q{r, n') , (5.9) 

where w(r; n, n') is the differential cross-section of the scattering n — > n' at the point 
r (for the isotropic scattering w(r; n, n') is equal to l/r(r)). 

MK conjectured further that the derived a-model makes sense also in the limit 
T —>■ oo, when it describes the clean system. Let us note, however, that one cannot 
simply set = {i.e. remove the disorder completely). Indeed, then we would 
get a particular clean system with uniquely defined energy levels, so that the DOS 
will be a sum of (5-functions. This is certainly not what we want (or what we are 
able) to calculate in the cr-model approach. The correlation or distribution functions 
that we are discussing imply necessarily some averaging. The MK derivation seems 
to remain meaningful if Ar ^ 1. In this case the disorder is strong enough from the 
quantum point of view, i.e. we study not a particular quantum system but rather a 
large ensemble of systems. On the other hand, the condition that the disorder does 
not influence the classical dynamics is r ^ L/vp, where L is the system size. Since 
L/vpA ~ {k-pLY^^ ^ 1 in the semiclassical limit, the double inequality L/v-p <^ 
can be satisfied. It means that the disorder is classically negligible, while it mixes 
strongly energy levels of the quantum system. 

^•^As has already benn mentioned, MK obtained the baUistic action in a different (Wess-Zumino- 
Witten) form, which is however equivalent to (^.61). 
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In a diffusive sample the density relaxation is governed by the diffusion operator 
with eigenvalues > (and eo = corresponding to the particle number conserva- 
tion in a closed sample). The positiveness of corresponds to the exponential decay 
of a density perturbation with time. What are their counterparts 7^ in a chaotic 
sample? One might naively think that, since the evolution operator e~^* is unitary, 
the corresponding eigenvalues e"'''''* are of absolute value unity, so that 7^ are purely 
imaginary. This is, however, incorrect. It is known that an ultraviolet regularization 
(projection onto the subspace of smooth functions) shifts all 7^ (except 70 = 0) from 
the imaginary axis, giving them a positive real part, Re7^ > 0. These eigenvalues are 



known as Ruelle resonances ||129|| . The corresponding regularized evolution operator 
(called Perron- Frobenius operator) describes irreversible classical dynamics in a chaotic 
system. The above ultraviolet regularization may be physically understood as an in- 
finitesimally weak noise introduced in the system to make the dynamics irreversible. 
The identification of the eigenvalues 7^ determining the non-universal corrections to 



the spectral statistics in a chaotic system with the Ruelle resonances was done in [|130| 



5.2.3 cr-model from energy averaging. 

In a subsequent paper, Andreev, Agam, Simons, and Altshuler (AASA) ||131|| proposed 
another derivation of the ballistic a-model. Instead of averaging over disorder, they 
started from a completely clean system and performed the energy averaging in a certain 
spectral window. After the Hubbard- Stratonovich transformation, they arrived at the 
action of the form 

S[Q] oc Strr{-2f -H[H, f]A - iujQA) , (5.10) 

where H is the quantum Hamiltonian, T and Q = TAT^^ are operators in the Hilbert 
space (and have on top of this the usual supermatrix structure), and Strr includes the 
supertrace over the supermatrix indices and the trace over the Hilbert space. The next 
(and crucial) step is the semiclassical expansion. Going to the Wigner representation 
and using the fact that in the semiclassical limit the Wigner transform of a commutator 



is a Poisson bracket of the corresponding Wigner transforms, AASA reduced ( |5.1CI| ) to 



the form (|5.6|) . More recently Zirnbauer demonstrated, however, that an additional 



averaging is necessary to guarantee the condition of a slow variation of the Wigner 
transform T(r, p) required by the semiclassical expansion. 

More specifically, Zirnbauer considered the level correlations for a unitary map 
(rather than for a Hermitian Hamiltonian). Similarly to the AASA approach, he aver- 
aged over the quasienergy e**^. This can be done via the "color-flavor transformation" 



19 



j d0expi(^le*'^^i + ^^6-^*^2) = J d/i(Z, Z) exp(^lZ^2 + ^2^^i) , (5-11) 
where ^'i, are supervectors, ^ ~ ^ 2 1 ] matrix from the coset space 
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and dfi{Z, Z) is the corresponding invariant measure.^ This transformation replaces 
(in the case of unitary maps) the energy averaging and the Hubbard-Stratonovich 
transformation. Moreover, the saddle-point approximation is not needed in this case, 
since the r.h.s. of ( |5.11| ) is aheady an integral over the coset space. Zirnbauer showed 



further that, in order to justify the semiclassical expansion, one has to average over an 
ensemble of maps, 

f/(0 = exp (^zJ2jkXk/hj U , (5.12) 

where C,k are random variables with the disorder strength scaling as {^D oc < a < 
1, in the limit ^ — ^> 0. Since a > 0, all the members of the introduced ensemble have 
the same classical limits. On the other hand, the condition a < 1 ensures that the 
disorder is strong from the quantum point of view. This procedure is thus physically 
similar to the derivation of MK with L/fp ^ r ^ A^^ (see above). 

5.2.4 Non-universal corrections and statistical noise. 

Let us discuss now the implications of the ballistic cr-model (|5.(j|) for the level statistics. 



Since the system is assumed to be chaotic, the only zero-mode [i.e. a field Q(r) yielding 
zero when substituted in the first (kinetic) term of ( p. 6] )] is Q = const, so that in 
the zero-mode approximation the OD cx-model and thus the RMT are reproduced, as 
expected.]^ Corrections to the RMT have the same form as discussed in Sec. 
p.4| , with the diffusion eigenvalues replaced by the Perrron-Frobenius eigenvalues 
7^. In a strongly chaotic system of a characteristic size L, a typical relaxation time 
(the counterpart of the diffusion time ^d) is of the order of the flight time, ~ L/vp. 
Therefore, the ballistic counterpart of the dimensionless conductance g ^ 1/ Ato can be 
estimated as ~ l/Ate ~ {kpLY^^ ~ jsjid'-'^)/^ ^ where ~ E//S. is the level number 
around which the statistics is studied. In particular, for the most often considered case 
of a 2D billiard g ~ A^^/^. We will discuss the deviations in more detail in Sec. 
where important differences compared to the diffusive case will be demonstrated. 



Prange [ |132|| has discussed conditions of observability of the non-universal correc- 
tions to the spectral form-factor K{t) around r = 27r predicted by the cr-model. He 
pointed out that K{t) is a strongly fluctuating function with the r.m.s. deviation equal 
to the mean value. After averaging over a window of width ~ 1/(7 around r = 27r the 
noise amplitude is reduced to {g/NY^"^. Therefore, to detect the deviation SK2n{T), 
whose amplitude is ~ 1/g, the following inequality should be satisfied: 

0x1/2 T 

^) «- =^ N:^g'. (5.13) 

However, as we have just discussed, g ~ N^^'^ for a generic 2D chaotic billiard, so 
that the condition ( p.l3| ) is not fulfilled. Therefore, the non-universal correction to 



and Z are identical to —W12/2 and —W21/2 respectively if the parametrization ( ^.22 ) is used 



34 j 

^^For an integrable system any function Q depending on integrals of motion only would be a zero 
mode, thus invalidating this consideration. 
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the spectral form-factor is not observable for an individual quantum system, when the 
only averaging available is the energy averaging. This points again to the necessity 
of the additional averaging over an ensemble of quantum systems having the same 
classical limit ||132|| . Let us note that for a diffusive system, Eq. ( [5.13[ ) can be satisfied 



without problems. In this case g and are two independent parameters, and one can 
consider the limit of arbitrarily large at fixed g. Therefore, it is in principle possible 
to extract the non-universal correction from a single disordered sample by using the 
energy averaging. 

5.2.5 Problem of repetitions. 

The following subtle point concerning the non- universal correction to R2{uj) predicted 
by the ballistic cr-model is worth mentioning here. The smooth (Altshuler-Shklovskii) 
part of the correction, Eqs. ( |3.31D , ( p.57[ ), can be written as follows 



/-oo 



27r2 



Re / dtte"^' Tie- , (5.14) 







which is easily checked by using the fact that the eigenvalues of the evolution operator 
e~'^* are equal to e"'*'''*. Expressing the trace in terms of the Gutzwiller sum over 
periodic orbits, one gets [[L33| , |134|| 

p r=l ' ^ P ' ' 

where the index p labels primitive orbits with periods Tp, the summation over r takes 
into account repetitions of the primitive orbits, and Mp is the monodromy matrix 
characterizing the dynamics in the vicinity of the orbit p. Substitution of ( |5.15| ) into 
( |5.14| ) readily yields 

41sH = VT,^f . (5.16) 

2.ASV ; 27r2 ^ p det M^I - 1 I ^ ^ 

p J.— I I \ p J 



On the other hand, one can calculate the two-level correlation function semiclassically, 
by starting from the Gutzwiller trace formula for the DOS ||135| , |134|| 



4^ = 1 + Re^ V T, V ^^P{J[^.(^) - ^W2]} 
iv) TT ^ det M^-l 1/2 ' ^ ' 

where Sp is the action and fip the Maslov index of the orbit p. The two-level correlation 
function is thus represented as a double sum J2pp' o^^r the periodic orbits. Assuming 
that the terms with p ^ p' vanish upon energy averaging due to randomly varying 



phase factors ("diagonal approximation" introduced by Berry), one finds |47 



d(c) A _ r>^ ^2 V" exp(^u;Tpr) 

^2,diag('^) - ^^^l^^v 1. |det(M;-l)| • ^^-^^^ 
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Comparing equations (|5.16| ) and ( p.l8| ), we see that they almost agree, the only dif- 



ference being in the extra factor r in ( |5.18|) . This discrepancy between the results 



of the ballistic a-model and the diagonal approximation was emphasized by Bogo- 
molny and Keating ||136|| . The same semiclassical analysis was performed earlier by 
Argaman, Imry, and Smilansky in the context of diffusive systems. In that case, 
however, the relevant periodic orbits determining the non-universal behavior have a 
length ~ vptj) ^ vpT, i.e. they are much longer than the shortest periodic orbits. For 
this reason and in view of the exponential proliferation of the primitive periodic orbits 
with length, one can neglect the repetitions, keeping only the r = 1 term in the trace 
formulas. Then Eq. ( |5.16|) and (|5.18|) are in full agreement with each other. On the 



other hand, in the ballistic case the shortest orbits of the length ~ L are relevant, and 
there is no parameter which would justify neglecting the repetitions. There exists thus 
a real discrepancy between the two formulas, which still awaits its explanation. 

5.3 Billiard with diffuse surface scattering 



According to Sec. f>.2[ the level and eigenfunction statistics of a clean chaotic system 
(with an ensemble averaging discussed above) are described by the formulas of Sec. ^ 
^ with the Perron- Frobenius operator substituted for the diffusion operator. However, 
straightforward application of these results to a given chaotic billiard is complicated by 
the fact that the eigenvalues of the Perron- Frobenius operator are unknown, while its 
eigenfunctions are extremely singular. For this reason the cr-model approach has so far 
failed to provide explicit results for any particular ballistic system. In this subsection, 
we consider a ballistic system with surface disorder (a rough boundary) leading to 
diffusive scattering of a particle in each collision with the boundary ||13?| , [L38|] . Since 



the particle loses memory of its direction of motion after a single collision, this model 
describes a limit of an "extremely chaotic" ballistic system, with typical relaxation time 
being of order of the flight time. (This should be contrasted with the case of a relatively 
slight distortion of an integrable billiard [p.39| , |140|| .) This is a natural problem to be 
studied by the ballistic cr-model approach. While the assumption of the diffuse surface 
scattering makes possible an explicit analytic treatment of the problem ||137| , |138|| P^ 



the obtained results seem to reflect generic features of ballistic systems in the regime of 
hard chaos. To simplify the calculations, a circular geometry of the billiard is assumed. 
A similar problem was studied numerically in Ref. ||142|| for a square geometry. As 



usual, we consider the case of unitary symmetry; generalization to the orthogonal case 
is straightforward. 

Inside the billiard, the motion is free and the Liouville operator C is given by 
Eq. (15.31). Clearly, it should be supplemented by a boundary condition, which depends 
on the form of the surface roughness. As a model approximation we consider purely 



diffuse scattering ||143| , |144|| for which the distribution function (p{r, n) of the outgoing 



Very recently, the same approach was used [ 141 1 to calculate the persistent current in a ring with 
diffusive scattering. 
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particles is constant and is fixed by flux conservation^: 

ip{r,n) = 7T (Nn)^(r,n)dn, (Nn) < 0. (5.19) 

i(Nn')>0 

Here the point r lies at the surface, and N is an outward normal to the surface. This 
boundary condition should be satisfied by the eigenf unctions of £. The boundary 
condition breaks the naive anti-hermiticity of C, and all its eigenvalues (except the 
zero one) acquire a positive real part, as expected for the Perron- Frobenius operator 
of a chaotic system. 

Specifically, the eigenvalues 7 of the operator C corresponding to the angular mo- 
mentum / obey the equation 

1 r 

J,(0 = -1 + -/ desmeexp[2iW + 2^sme] = 0, (5.20) 



where ^ = R'j/vf, and R is the radius of the circle. For each value of / = 0, ±1, ±2, . . . 
Eq. (|5.20| ) has a set of solutions ^ik with = ^-z,fe = ^f„fc, which can be labeled with 
A; = 0, ±1, ±2, . . . (even I) or k = ±1/2, ±3/2, . . . (odd /). For / = = we have 
^00 = 0, corresponding to the zero mode (f{r, n) = const. All other eigenvalues have a 
positive real part IieC,ik > and govern the relaxation of the corresponding classical 
system to the homogeneous distribution in the phase space. The asymptotic form of the 
solutions of Eq. ( p.2(]| ) for large \k\ and/or \l\ can be obtained by using the saddle-point 
method, 

_ r 0.66/ ± 0.14 In/ ± 0.557rifc, < < / 
^^'"^ \ {\nk)/4 + 7rt{k + l/8), < / < A; " ^' ^ 

Note that for k = all eigenvalues are real, while for high values of k they lie close to 
the imaginary axis and do not depend on / (see Fig. |l|). 

5.3.1 Level statistics. 

The characteristic frequency separating the regions of the close-to-RMT and the fully 
non-universal behavior (analog of the Thouless energy) for the considered problem is 
u ~ tg^ ~ vp/R. In the low-frequency range, u Vp/R, the level correlation function 
is given by Eq. ( |3.48|) , where we expect A ~ 1/(7^ ~ (RA/vp)^. Indeed, the calculation 
yields p8| 

(RA \ ^ 
, a = -19/27 - 1757rVll52± 64/(97r2) ^ -1.48. (5.22) 
TTVp J 



^^The exact form of the boundary condition depends on the underlying microscopic model. In 
particular, the diffuse scattering can be modelled by surrounding the cavity by a disordered layer with 



a bulk mean free path I and a thickness I. The corresponding boundary condition |145, lAt \ differs 



from Eq. ( 5.19| ) by a parameterless function of order unity. For a review of the boundary conditions 



corresponding to various microscopic realizations of the rough surface see 1 147 
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40.0 




Figure 1: First 11 x 11 (0 < /c, / < 11) eigenvalues of the Liouville operator L in units 
oivplR-i as given by Eq. (|5.20|) . From Ref. ||137 |. 
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In contrast to the diffusive case, this constant is negative: the level repulsion is en- 
hanced with respect to the RMT. 

The high-frequency behavior of R2{'^) is expressed in terms of the spectral deter- 
minant ( |3.56| ) by the formulas of Sec. |3.3.3| . Calculating the spectral determinant for 



the present problem, we find |137| , |13^ 

D{s) = s-^ J] {l-isAhkiY\l + isAhki)-^ 

(2) N n Ji{isN-y^)Ji{-isN-y^) ' ^^'^^^ 

where N = {vp/RAY = {ppR/2)'^ is the number of electrons below the Fermi level. 
For high frequencies uj ^ vp/R this yields the following expression for the smooth 
(Altshuler-Shklovskii) and the oscillating part of the level correlation function: 

„ , , (ARVf Vf \V2 f ojR Ti\ 



It is remarkable that the amplitude of the oscillating part ( |5.25|) does not depend on 



frequency. This is in contrast to the diffusive case, where in the high-frequency regime 
{uj above the Thouless energy) the oscillating part i?osc(^) is exponentially small, see 
Eq. (|33|). 

5.3.2 The level number variance. 

The smooth part of the level correlation function can be best illustrated by plotting 
the level number variance ^2(5) [see Sec. p.4.2|| . A calculation [|137|| gives for s <^ N^^"^ 

7r2E2(s) = 1 + C + ln(27rs) + as'^/{2N) (5.26) 

and for s > iV^/^ 

.^E2(.)^l + C + ln-^--(--j cos(^--j. (5.27) 

Here C ~ 0.577 is Euler's constant, and the numerical constant a is defined by 
Eq. ( |5.22| ). The first three terms at the r.h.s. of Eq. ( [5. 26] ) represent the RMT contri- 



bution (curve 1 in Fig. |^). As seen from Fig. |^, the two asymptotics (|5.26 ) and ( ^.27] ) 



perfectly match in the intermediate regime, s ~ N^/'^. Taken together, they provide 
a complete description of ^2(5). According to Eq. (|5.27|) , the level number variance 
saturates at the value Sg^"* = 7r^^(l + C + ln(16A^^''^/7r^)), in contrast to the behavior 
found for diffusive systems [see Eq. ( p. 701 )] or ballistic systems with weak bulk disor- 
der [p.48| , |149|| . The saturation occurs at energies s ~ or in dimensionful units 
E ~ vp/R. This saturation of ^2(3), as well as its oscillations on the scale set by 



short periodic orbits, is expected for a generic chaotic billiard It is also in good 
agreement with the results for ^2(5) found numerically for a tight-binding model with 
moderately strong disorder on boundary sites [|142| . 
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Figure 2: Level number variance T,2{E) as a function of energy; s = E/A. Curve 1 
shows the RMT result, while curves 2 and 3 correspond to asymptotic regimes of low 
( |5.26| ) and high ( [5.27| ) frequencies. The saturation value Eg^"* is given in the text. From 
Ref. [T^. 
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5.3.3 Eigenfunction statistics. 

Let us recall that the system-specific contributions to fluctuations and correlations 
of the eigenfunction amplitudes were described in the diffusive case in terms of the 



diffusion propagator n(r,r'), see Sec. [4.5.3| . Its ballistic counterpart IIb is deflned 
as follows 



nB(ri, r2) = J dnidns g{ri, ni; ra, ng), 

Cgir,, m; r^, n^) = {nu)-' [5{r, - r2)6{n, - 112) - V''] . (5.28) 
Equivalently, the function nB(ri,r2) can be defined as 

POO P 

nB(ri,r2) = y dty dni^(ri,ni,t;r2) , (5.29) 

where g is determined by the evolution equation 

d 



+t;FniViJ ^(ri,ni,t;r2) = , t>0 (5.30) 

with the boundary condition 

g\t=o = {niy)-'[6{r, - ra) - V-']. (5.31) 



Equation ( |5.28 ) is a natural "ballistic" counterpart of Eq.( 3.42| ). Note, however, that 
nB(ri, i'2) contains a contribution ng^(ri, r2) of the straight line motion from r2 to ri 
(equal to l/(7rpF|ri — r2|) in 2D and to 1/2(pf|i'i — i'2|)^ in 3D), which is nothing but the 
smoothed version of the function /|(|ri — r2|). For this reason, n(ri, r2) in the formulas 
of Sec. ^ should be replaced in the ballistic case by n(ri, r2) = nB(ri, r2) — (ri, r2). 
At large distances |ri — r2| ^ \f the (smoothed) correlation function takes in the 
leading approximation the form 

V^mr,)\'mr,)\' = 1 + |nB(ri,r2). (5.32) 



For related results obtained in the semiclassical approach see Refs. ||150| , |151| 



Equation ( p.32| ) shows that correlations in eigenfunction amplitudes in remote 
points are determined by the classical dynamics in the system. It is closely related 
to the phenomenon of scarring of eigenf unctions by the classical orbits [[152| , |153|| . In- 



deed, if ri and r2 belong to a short periodic orbit, the function nB(ri, r2) is positive, so 
that the amplitudes |'?/'fc(ri)p and \ipk{T^2)\'^ are positively correlated. This is a reflection 
of the "scars" associated with this periodic orbit and a quantitative characterization of 
their strength in the coordinate space. Note that this effect gets smaller with increasing 
energy E of eigenfunctions. Indeed, for a strongly chaotic system and for |ri — r2| ~ L 
[L being the system size), we have in the 2D case nB(ri,r2) ~ )^f/L, so that the 
magnitude of the correlations decreases as 
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For the model of a circular billiard with diffuse surface scattering a direct calculation 
gives ||137| 



nB(ri,r2) = ni(ri,r2)+n2(ri,r2), (5.33) 

ni(ri,r2) = ni?)(ri-r2)-l^-^y"dVi<)(r'i-r2) (5.34) 
-V~' I dh',U^^\r, - r'2) + I dh\dh',U^^\r\ - r',); 

where n3^(r) = l/(7rpi7'|r|), and {r,6) are the polar coordinates. This formula has 
a clear interpretation. The function IIb can be represented as a sum over all paths 
leading from ri to r2, with possible surface scattering in between. In particular, IIi 
corresponds to direct trajectories from ri to r2 with no reflection from the surface, while 
the contribution 112 is due to the surface scattering. The first term in the numerator 
Ak"^ — 1 comes from trajectories with only one surface reflection, while the second sums 
up contributions from multiple reflections. 



5.4 cr-model for the kicked rotor. 

Up to now we have considered autonomous systems, with the Hamiltonian H having no 
time dependence. The supersymmetric cr-model approach is, however, also applicable 
to periodically driven systems {H{t) periodic in time t), as we are going to discuss. The 
standard system of the latter type is the quantum kicked rotor (QKR) |p,54| , |155| , pO| 
defined by the Hamiltonian 

H = — + ~kcose J2 ^it-niT). (5.35) 

m=— 00 

Here / = —ihd/dO is the angular momentum operator conjugate to the angle variable 
9] I is the moment of inertia, T the kick period, and k the kick strength. The classical 
version of this problem is characterized by a single dimensionless parameter K = kr/I. 
For a sufficiently large K the motion becomes globally chaotic and the system exhibits 
unbounded diffusion in the phase space (in the direction of the angular momentum), 
the diffusion coefficient being 



2t 



h2 

~ — for > 1. (5.36) 

AT 



In contrast to the classical problem, the quantum system depends non-trivially on 
both parameters k and T, since there are two dimensionless combinations which can 
be formed: k = k/h and r = hT/I. The classical limit corresponds thus to A; ^ 00, 
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r — s> at fixed kr = K = const. After tliis sliort reminder of the classical-quantum 
correspondence for the kicked rotor problem, we set ^ = 1 as in the other parts of this 
article. 

Numerical simulations of the QKR have shown that at long times the classical dif- 
fusion is suppressed, the phenomenon known as "dynamical localization". In Ref. | |156 | 



an analogy between this problem and the Anderson localization in ID systems was 
drawn. Later studies revealed a close connection between the QKR problem and the 
RBM ensemble (see Sec. 2.5.4 ). The evolution (Floquet) operator of the QKR has in 
the angular momentum representation the form 



exp\^-tU^^Ji_i,{k) , (5.37) 



where Jm{k) is the Bessel function. Since Jm{k) — (27rm)~^/^(efc/2m)"* for m ^ 1,/c, 
we see that the matrix elements of U are exponentially small for | / — /'| ^ k. Therefore 
the matrix U is indeed banded, with the bandwidth 6 ~ ^ 1. Furthermore, the 
second factor on the r.h.s. of ( ^.37] ) mimics a generator of pseudorandom numbers (of 
absolute value unity). Indeed, numerical simulations |3^, |1571| have demonstrated 
that the statistical properties of the two models (QKR and RBM) are very close to 
each other. 

On the other hand, the RBM ensemble was shown to belong to the class of quasi-lD 
disordered systems described by the ID a-model ||32[, see also Sec. ^4.3.1| . The solution 
of the ID cr-model made possible a detailed analytical study of eigenfunctions of such 
systems []7D[, see Sec. ^7^ . While the pseudo-RBM properties of the QKR evolution 
operator f/„„', as well as the numerical simulations, suggested strongly that these re- 
sults are also applicable to QKR, a formal derivation of this fact has been missing 
until recently. This gap was filled by Altland and Zirnbauer (AZ) |[158|| who achieved 
a mapping of the QKR onto the ID a-model. The idea of their calculation is essen- 
tially the same as the one described in Sec. ^.2.3| . Performing the averaging over the 
quasienergy spectrum with making use of the transformation ( |5.11| ) and carrying out 
the semiclassical expansion, AZ obtained the ID version of the cr-model ( p.l6| ) (or its 
orthogonal-class counterpart for unbroken time- reversal symmetry), with D equal to 
the classical diffusion constant of the rotor ( p.36| ) and u = T/2tt. This is precisely 
what could have been expected from the analogy with disordered wires or the RBM 
ensemble; the above value of u being the density of quasienergies ujk corresponding 
to the eigenvalues e*'^'=^ of the Floquet operator. This mapping allowed AZ to con- 
clude that the statistical properties of the QKR are identical to those of the quasi-lD 
disordered systems. In particular, the localization length Lioc of the QKR governing 
the decay of a typical eigenf unction, |'?/'^(0|typ oc exp(— |/ — /o|/-^ioc), is found to be 
Lioc = PttuD = {(3/8)k'^, in agreement with numerical simulations of Shepelyansky 

■^^Shepelyansky considered also another localization length - that of a steady state - and concluded 
that it is larger by a factor of 2. This is in disagreement with the field-theoretical results, which give 
the same value Lioc = (/3/8)fc^ for the both lengths. Presumably, the discrepancy is due to insufficient 
numerical accuracy of evaluation of the steady state asymptotic decay rate for ^ 1 in |159|. 
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It should be noted, however, that the problem of sufficiency of the energy averaging 
discussed in Sec. |5.2.3| is equally apphcable here. In ||160|] AZ acknowledged that an ad- 



ditional averaging over an ensemble of rotors having the same classical limit {i.e. of the 
type described in Sec. [5.2.3| ) was needed to justify the semiclassical expansion which is 



in the heart of the derivation in [|158[|. In fact, the conclusion of the existence of such an 



implicit ensemble averaging in ||158|| can also be supported by the following argument. 
Following ||158|| , one can calculate, e.g., the distribution function of the quantity of 
the type ( [4.31| ), v oc \'ipai^i)'ipai^2)\, which will have the log-normal distribution (|4.32|) . 



In particular, the far tail of this distribution is crucially important for identifying the 
localization length of a typical eigenstate on the basis of the average value (v) (which 
has a 4 times smaller decay rate). However, to be able to find the whole distribu- 
tion, one should average over an exponentially large exp(|/i — /2|/4Lioc)] number 
of eigenfunctions, while the energy averaging alone reduces effectively to an averaging 
over ~ \h — hi eigenfunctions only and is thus insufficient. 

5.5 Concluding remarks. 

We finish these notes by comparing briefly the supersymmetric and the semiclassic 
(periodic orbit) approach to the spectral statistics of chaotic systems. The main ad- 
vantage of the supersymmetry method is that it allows one to get the RMT results in 
the leading approximation. In contrast, the semiclassical approach is only justified for 
times much shorter than the Heisenberg one, r/27r ^ 1. Obtaining the GUE results 



within this approach requires using an ad hoc regularization prescription ||136| , |161 
Even the problem of calculating the perturbative (in r) corrections to the leading be- 
havior K{t) = r/vr in GOE (for GUE the perturbative corrections are identically zero) 
has not been solved semiclassically. Therefore, the cx-model approach is the only known 
method allowing to obtain the RMT results in a controlled way. 

On the other hand, the ballistic cr-model approach is also not free from problems. In 
particular, construction of a regular expansion inl/g has not been achieved in this case 
(in contrast to the diffusive cr-model). While the leading non-universal contributions 
(discussed above) to the energy level and eigenfunction statistics are only determined 
by the density relaxation modes 0^ and their eigenvalues 7^, the higher-order (in l/g) 
terms are induced by interaction of these modes. In the case of the diffusive cr-model, 
the corresponding contributions are known as weak localization corrections and can 
be calculated systematically. At the same time, recent attempts to perform such a 



calculation in the ballistic case ||162| , |163|| did not produce unambiguous results, because 
of the singular nature of the expressions obtained. This problem (as well as that of 
repetitions, see Sec. [5.2.5| ) remains a challenge for future research. 
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